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ABSTRACT 


The application of Direct Bit Detection methods are analyzed and evaluated in 
digital communication systems employing Differential Phase Shift Keyed (DPSK) 
modulation. Assuming an additive white Gaussian noise interference model, and direct 
bit detection methods used in the receiver, 8-DPSK and 16-DPSK communication 
systems are considered and their performances evaluated in terms of the delivered bit 
error rate as a function of signal-to-noise ratio. The advantages and disadvantages of 
Direct Bit Detection Receiver (DBDR) systems used in conjuction with differential 
phase encoding are determined, with specific application to 8-DPSK and 16-DPSK 
modulated signals. The effect of improper receiver carrier phase synchronization 1s 
considered in detail, and resulting performance degradations are evaluated. Numerical 
results show that for receiver phase errors of more than a few degrees, severe 
performance degradations result for both DPSK modulation schemes, unless a 


complete phase reversal (i.e.,180°) takes place. 
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I. INTRODUCTION 


M-ary differential phase shift keved (M-DPSK) modulation is a widely used 
technique in digital communication applications. Bandwidth efficiency, relatively good 
noise immunity, constant signal envelope, and simplicity of implementations, make this 
scheme particularly attractive for use over satellite, terrestrial radio and voiceband 
telephone channels. While system analyses pertaining to the performance of M-DPSK 
receivers abound in the literature, treatment is usually restricted to the case where 
signal transmission takes place over an additive white Gaussian noise channel. [Ref. 1] 

Since the digital information transmission is based on carriers modulated by the 
symbol waveforms, the symbol error probabilities which quantify receiver performance 
can often be computed directly as demonstrated by the many derivations of such 
probabilities found in the literature ( see [Ref. 2: pp. 228-234] and [Ref. 3: pp. 200-212] 
for example ). However, from the point of view of the recipient of the digital 
information, it is the Bit Error Ratio (BER) that becomes important in digital 
communication applications. While results on the BER performance of M-arv 
modulation receivers is far more limited, closed form expressions for the BER of M- 
PSK receivers have been obtained by Lee [Ref. 4], with similar results independently 
derived by Tan [Ref. 5]. 

This thesis is devoted to analyzing the application of direct bit detection methods 
which allow the recovery of bits in the received symbol individually, in differential 
PSK digital communication systems. For analysis purposes, 8-DPSK and 16-DPSK 
svstems are considered and their performances evaluated. The particular question to 
be answered is whether for 8-DPSK and 16-DPSK modulated signals, direct bit 
detection methods provide comparable performance to conventional phase detection 
receivers, especially under conditions of improper receiver carrier phase 
synchronization. 

In this introductory chapter, background information on M-PSK communication 
systems is given. Known optimum recelver structures are presented, and Direct Bit 
Detection Receivers (DBDR) are introduced. In Chapters II and III a presentation of 
the transmitter and receiver structures respectively is given along with their operational 
properties, and the logic implemented in the DBDR for 8-DPSK and 16-DPSK 


I] 


modulated signals is discussed in detail. Chapter IV includes the noise performance 
analyses of 8-PSK and 8-DPSK receivers, and Chapter V presents similar receiver 
analyses for 16-PSK and 16-DPSK. In both of these chapters, receiver local oscillator 
phase errors have been incorporated so as to be able to account for degrading effects 
due to improper carrier phase synchronization. The obtained mathematical expressions 
of receiver performance are evaluated on a computer and graphically presented as 


plots of bit error ratio as a function of signal to noise ratio (SNR) and phase errors. 


A. BACKGROUND 

Direct bit detection methods utilize two signal processing channels in the receiver 
in order to recover the digital information without phase angle measurements. 
Moreover, these methods exhibit optimum performance in a BER sense when carrier 
svnchro..ization is completely achieved. Because of the attractiveness and simplicity of 
DBDRs, some studies have been carried out in this area, and a modified two-channel 
receiver for 8-PSK has been built and its noise performance measured by Thompson 
(Ref. 6]. Furthermore, noise performance evaluations of DBDR for 8-PSK have been 
presented by Myers and Bukofzer [Ref. 7], while similar analyses for the 16-PSK case 
can be found in Reference 8. Advantages of direct bit detection methods are well 
documented in these referenced works, in which all analyses have been done assuming 
no receiver local oscillator phase error, Gray Code bit to symbol mapping and an 
additive white Gaussian noise interference model. We present here 8-PSK, 8-DPSK, 
16-PSK and 16-DPSK communication systems under similar assumptions, except now 


the receiver has an assumed Known local oscillator phase error. 


B. THEORY OF M-PSK COMMUNICATION SYSTEMS 
1. Representation of M-PSK Signals 
For M-PSK modulation with equal signal energies, a convenient 


representation [Refs. 3: p.p. 192-205] of the signal set is given by 


s(t) = J2E,/T, cos(2mftt+ 2n(i-1)/M), 1=1,2,.,.M , 0 St ST, (1.1) 
bie 
E.=f s(t) dt , i= 1,2,..,8 (1.2) 
0 
T, = symbol duration 
fy) = carrier frequency (Hz.) 


for convenience > > alee is assumed. 


A suitable orthonormal signal set for the representation of the above signals is 


O (N= V2/T, cos(27f,t) 

9(t)= /2/T, sin(2zf,t) tS ies 63) 
so that by using trigonometric identities, s. (t) can be expanded in terms of @,(t) and 
@.,(t). That is, 


s(t) = @,(t) /2E, cos[ 2m(i-1)/M ]- @ ,(t) /2E, sin[ 2m(i-1)/M } (1.4) 


Where 1=1,2...... VM. 
A plot of the signal constellation for M-PSK is shown in Figure 1.1 for 


th 


various values of M, where the coordinates of the i“ signal vector representing the 


signal s.(t) are 
{ /2E, cos[2m(i-1)/M]} ,- { /2E, sin[2n(i-1)/M] } (1.5) 


For M-DPSK, the signal representations are exactly same as those already presented, 
except that the digital information is transmitted in terms of phase differences between 
consecutive signal transmissions, rather than absolute signal phases. 
2. Gray Coding for M-PSK Signals 
Before proceeding with the Gray coding for M-PSK signals, we introduce 
some definitions which will henceforth be used 


1. Symbol : Transmitted or received signal in sinusoidal form (s(t) where 
i= 1,2,...M) which is related to the symbol state. 


2. Symbol state : The binary digit assignment to the symbols which always 
consists of a block of log, M bits in length(S, where i= 1,2,...,M.) 


3. Data state : Block of log,M bits corresponding to the data that must be 
delivered to the intended user, and which is processed at the encoder input to 
determine the symbol to be transmitted. (D. where i= 1,2,....,.M.) 


The relationship between symbol error probability and bit error probability in 
the corresponding log,M bit groups depends upon the assignment of bits to the symbol 
states. The preferred bit to symbol state mapping in most cases is the Gray Code, in 
which the symbol states which correspond to adjacent phases differ in only one bit 


position. This is because in the coherent demodulation process the most likely error 
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Figure 1.1 Signal Constellations for (a) BPSK (b) QPSK (M=4) (c) 8-PSK (d) 16-PSK. 


involves mistaking a correct symbol with its nearest neighbor and so such symbol 
errors translate into single bit errors. 

In Direct Bit Detection Receivers (DBDR), Gray code assignment to the 
symboi states are selected for the above given reason. In this thesis, a modified Gray 


code is used, which for comparison purposes is shown along with the conventional 


TAB | 


COMPARISON OF CONVENTIONAL AND MODIFIED GRAY 
CODING SCHEMES 


States. Coding Coding 
8-PSK 16-PSK $-PSK 16-PSK 
co 000 0000 HI M11 
Ss 001 0001 101 1101 
Ss O11 0011 001 1001 
Sh 010 0010 O11 1011 
Se 110 0110 010 0011 
oF M1 O11] 000 0001 
Ss 101 0101 100 0101 
Ss 100 0100 110 O11 
s 1100 0110 
oe 1101 0100 
S HH 0000 
o 1110 0010 
a 1010 1010 
oh, 1011 1000 
Ey. 1001 1100 
oy 1000 1110 


Gray Code in Table 1. 
Using the modified Gray Code on Table 1, the vectors representing the 
symbols with components along the @, and @, axes have been shown (see Figure ZZ) 


such that r, is the symbol component along the 9, axis and r, is the symbol 
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component along the @, axis. Therefore for M=8, (see Figure 1.2.a) in order to 
recover Most Significant Bit (LB) of the transmitted symbols correctly, r, >0 whenever 
the symbol states §,, S85, S, and §, are transmitted and r, <0 whenever the symbol 
states §;, S,, S; and §, are transmitted. The Middle Bit (MB) of the transmitted 
symbols can be recovered correctly if Cre - (5)? > 0 whenever the symbol states §,, 
S,, 8, and S. are transmitted and ce - (as < 0 whenever the symbol states §., S,, 
S, and §. are transmitted. The Least Significant Bit (RB) is recovered correctly if 
r,>0 whenever the symbol states $,, S,, $3, and §, are transmitted and r,<0 
whenever the symbol states $<, 8,, $5, and S, are transmitted. 

Similar logic has been set up for the M=16 case (see Figure 1.2.b) but now 
each symbol is related to 4 bits in a symbol state. Therefore for this case, 4 bits must 
be recovered from each received svw™bol, where the bits in a symbol state are labeled 
Most Significant Bit (LB), Middle Bit Right (MBR), Middle Bit Left (MBL), and Least 
Significant Bit (RB). For example from Figure 1.2.b, §. is O101, corresponding to 
LB=0, MBL=1, MBR=0 and RB=1. In order to recover the LB correctly, it is 
necessary that r,>0 for transmitted symbol states §,, $5, 83, 84, 813, 8,4, 3,5, and 
Sj¢ and that r, < 0 for transmitted symbol states §., S,, S., 8, S89, 819, 81), and §)5. 
Te r,> QO whenever the symbol states Si) S>, S33 34; Ss, Ses S5, and S. are transmitted 
and r, <0 whenever the symbol states So» Sigs Siy> Sia» S13» S14 S159 and Sig are 
transmitted, the RB will be recovered correctly. If CG) for transmission of 
symbol states §,, 85, 85, 83, So, Sig, 8,5, and S)¢, and (= ye for transmission of 
symbol states Do, arias, ae Siy> Sia» S13 and S14 the MBL will be recovered 
correctly. Finally, if [ ee - (iA ee se > 0 upon transmission of the symbol 
states $1, Sq > Scr Sqy Sq: Sip) Sign andes ap 2 yr) ne 
transmission of the symbol states on the states 85, 83, 8, 85, Sig, S11» Sy 4 and 815, 
then the MBR will be recovered correctly. 

ln subsequent chapters, the receiver logic and its performance will be 
explained in detail. Additionally, application of this methodology will be presented in 
so far as differential PSK modulation is concerned. Before doing so however, a 
summary of Known results on the performance of conventional and DBDR’s for M- 
PSK and M-DPSK is presented next. 

3. Known Results on the Performance of M-PSK and M-DPSK Receivers 

In this section, known results on the receiver performance for M-PSK and 

DPSK modulated signals are presented for two cases, namely conventional receivers 


and Direct Bit Detection Receivers (DBDR). 
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Figure 1.2 Symbol State Assignments for (a) 8-PSK (b) 16-PSK. 


a. Conventional M-PSK and M-DPSK Receivers 

Coherent demodulation of M-PSK implies in principle the use of M signal 
processing channels in the receiver, as shown in Figure 1.3 (from Ha [Ref. 10] ) where 
different forms of the optimum receiver structure of M-PSK signaling have been given. 
Performance analysis results for optimum receivers are well-documented in terms of 
symbol error rate as a function of signal to noise ratio [Refs. 2,9], while more recent 
analvses have focused on the receiver probability of bit error (see [Ref. 4] and 
[Ref. 3: p.p. 204-207 |] for example). Svmbol error probabilities can be computed from 
knowledge of the channel and signal characteristics. However for transmission of 
binary data, when analyzing system performances as a function of different levels of 
modulation, the bit error probability rather than the symbol error probability is of 
interest, as previously explained. In Prabhu [Ref. |: p.198], the bit error prebutulity 
has been presented only for the case of M-ary orthogonal signal sets. However, since 
the M-PSK signal set is not orthogonal, those results are not applicable here. The 
closed form expression for the BER of M-PSK has been determined, when a Gray 


Code bit mapping is used [Ref.2: p.198]. The results can be summanized as follows 


n 
P.M) = » (k/n) P, (1.6) 
k=] 
where 
M = Number of distinct signal waveform, 


P, = Probability of k bits in error in a received n bit data 
block [Ref. 5: pp.24-32], 

n = log, M, 

P.(M) = Bit error probability. 


A derivation of Equation 1.6 and its numerical evaluation (see Figure 1.4 ) has been 
carried out by Tan [Ref. 5: p.33]. For M-PSK signalling, receiver performance analyses 
in terms of bit error rate have also been presented by Lee [Ref. 4: p.491], and plotted in 
Fae uvert 25. 
b. Direct Bit Detection Receivers (DBDR’s) 

DBDR structures for M-PSK have been derived for various values of M 
[Refs. 6,7,8]. As previously indicated, this method is particularly well suited for M-PSK 
communication systems. The DBDR structures for 8-PSK (see Figure 1.6) and 16-PSK 
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Figure 1.3. Block Diagram of Optimum Receiver (a) Correlator Realization 
(b) Matched Filter Realization (from Ha (Ref. 10]). 
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Figure 1.4 Bit Error Probabilities Versus E,/No for 16-PSK, 8-PSK and QPSK [Ref. 5 p.31]. 
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Figure 1.5 Probability of Bit error Versus E, {Ny [Ref.4]. 
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(see Figure 1.8) have been presented and their performances have been derived in 
References 6,7,8. The BER is quantified by the probabilities of a bit error, denoted PE. 
For 8-PSK, PE has been derived by Myers and Bukofzer [Ref. 7: p.13] 


where 


PE= 1/3 { Q2 YSNR cosa] + Q[2/SNR sina] + 
+ Q2 /SNR (cosa@+ sind)] + Q[ 2 YSNR ( cos@- sing) | - 


-2Q[2 YSNR (cosa - sing)]. Qf 2. “SNR (cosa+ sina) | } (1.7) 
where 


SNRe= E/2N j= 3ER NG 
E.= Symbol energy 
Ng = Noise power spectral density level (one sided) 


a = Basic signal phase angle (see Figure 1.2). 
and 


CO 


Q(x) = f [ 1//2m] exp(-u*/2) du (1.8) 


X 


The derivation of Equation 1.7 assumes an AWGN channel, Gray coded symbols and 
equal a priori bit probabilities (i.e. , 1/2). Receiver performance in BER sense is 
plotted in Figure 1.7 for various values of @ as a function of SNR. 

For the 16-PSK case, probability of bit error is derived by Bukofzer 
(Ref. 8] and the resulting receiver performance is plotted in Figure 1.9 for various 
values of the angles a and }, where B is shown on Figure 1.2.b. 

In both analyses, it has been assumed that the receiver's local oscillator 
does not have any phase error. In this research, similar analyses to those conducted 
in References 7 and 8 have been carried out by considering the effect of a receiver 
local oscillator phase error for 8-PSK and 16-PSK cases, as well as the noise 
performance of 8-DPSK and 16-DPSK DBDR, presented as BER as a function of 
SNR 
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Figure 1.6 Block Diagram for 8-PSK DBDR 
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ll. TRANSMITTER 


The generation of differenual PSK signals can be performed by pre-processing 
the data bits. This chapter is devoted to analyzing the M-DPSK transmitter structure. 
One general form of the transmitter in block diagram form for DPSK communication 
systems is shown in Figure 2.1. Explanation of the general elements in Figure 2.1 will 
be included along with the derivation of the transmitter structure for 8-DPSK and 
16-DPSK modulated signals. 


A/D 
Converter 


Encoder Carrier 


Modulation 





Figure 2.1 DPSK Transmitter Block Diagram. 


A. TRANSMITTER STRUCTURE FOR 8-DPSK AND 16-DPSK MODULATED 
SIGNALS 


The main task in generating the DPSK transmitter involves the design of the 
encoder, which depends on the algorithm used to relate data states to the symbol states 
prior to transmission as a sinusoidal waveform. Basically, the transmitted symbol is 
the quantized equivalent of the angle between the successively transmitted svmbols. 
For reasons previously explained, a Gray coding scheme is utilized in the svmbol state 
to symbol mapping for the two special cases of M-DPSK signalling considered here, 
namely M=8 and M=16. 


a 


The symbol differential phase relative to the data state to be transmitted is given 
on Table 2 for 8-DPSK and on Table 3 for 16-DPSK modulated signals. [In both 
Tables, it is assumed that an equiphase arrangement for the transmitted symbols is 


used. 


TABLE 2 


GRAY CODE REPRESENTATION OF PHASES FOR 8-DPSK SIGNAL 
Sie 


Differential 
Phase 


0 
@ 


Sy 
oo 
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As mentioned earlier, the main purpose of the encoder is to generate a mapping 
algorithm whereby the phase difference of consecutive symbol outputs is related to the 
input data in accordance with the information in Table 2 and Table 3. In order to 
achieve successive signal transmissions the following algorithm is implemented in the 
transmitter 


1. Store the present data state (D., i=1,2...,M) and the previous symbol state 
(g.,1= 1,2,3..,.M). 


=}? 
2. From the Table 2 or Table 3, find the corresponding differential phase value 


associated with the present data state as Na, where N is an integer with values 
of OF NSEM-I. 


3. Generate a symbol state, such that the corresponding symbol will have a phase 
angle Na@ degrees larger ( in additive sense, modulo 27) than the previously 
transmitted symbol. 
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For example, assuming the previous symbol state is S. (0101) and present data state is 
D, (9110) which corresponds to the differential phase value of 4a@ (see Table 3.). From 
Figure 1.2.b, proceed 4a degrres counterclockwise from §., resulting in the symbol 
state is $,4 (1111) will be the symbol output from the encoder. 

This algorithm can be shown in matrix structure, in which the row entries and 
the column entries of the matrix are the symbol states and data states respectively, and 
the elements of the matrix are the symbol state which correspond to the appropriate 


phase difference in the transmitted symbol (see Figure 2.3 and Figure 2.2). For 
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example, under the same assumptions given in the previous paragraph, by using Dy 
(0110) as a row entry and §. as a column entry in the matrix in Figure 2.3, present 
svmbol output from encoder is found as (1111) at the intersection of the row and 


column of the matrix. 


a 


SYMBOL STATE 


Sp Soo oS Smo 

DATA 000 100 110 111 101 001 O11 O10 
STATE 

Do= 000 000 100 110 11 10k CGO ore 

D.= 100 100 110 111 101 001 O11 010 O00 

D,= 110 110° 11) LOL 001 Ol .W10 1000 100 


Dee EL Dat LO O01 Ol Laer Ceo etree ro 
D,= 101 LOL 0017 Ol COLO OOO rfOCr reer 


D> O01 COOL OLIN OLOr OC MOC mire erin ot 
Dd, Or 011 010 000 100 110 111 101 O01 
D.= 010 010 000 100 110 111 101 001 O11 


Figure 2.2. Encoder Matrix for 8-DPSK Signals. 


The hardware implementation of the encoder structure can be accomplished in 
several forms. A practical approach to building the encoder can be achieved by the use 
of EPROMIs (Erasable Programable Read Only Memorv). Once the truthe table the 
programming of the EPROM is a trivial task to accomplish. Demultiplexers can also 
be used to build the encoder circuitry; however, because such a device would be hard 
wired, it would not have the reprogrammability of the EPROM implementation. 

In this research, the 8-DPSK and 16-PSK communication svstems are evaluated 
using a computer. Appendix A includes the source file which is used to generate the 
possible outputs of an M-PSK and an M-DPSK communication system. 

The carrier modulator uses the encoder outputs (symbol states) as an input, and 
generates the symbols (signals) suitable for transmission. The signal transmission 


involves sinusoidal voltages whose phase is related to the symbol states (encoder 
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outputs) and of carrier frequency suitable for the transmission channel. The signal 
generation can be achieved via a multiplexer, whose structure is shown in Figure 2.4. 
For this kind of implementation, (M-1) phase shift elements are needed with set phase 
shift values. Another approach would involve an implementation using a D/A 
converter along with a phase modulator which can be obtained with a VCO (see 


[Ref. 11: pp.269-271]. 


A Ae cos[2nf tt 


MULTIPLEXER 





Figure 2.4 Carrer Modulator Structure for M-DPSK Signaling. 


B. COMPUTER GENERATION OF THE TRANSMITTER OUTPUTS 

For computer evaluation of the M-arvy DPSK communication svstem , a 
FORTRAN program called BELIZ has been written (see Appendix A). The program 
generates all possible data state and symbol state combinations for M-DPSK, and 
produces the corresponding symbol states outputs. The program outputs are given in 
Appendix B in tabular form along with the transmitter matrix of Figure 2.3 and Figure 


2.2, specified in decimal form. 


Pit by EX 


A. RECEIVER STRUCTURE FOR 8-DPSK AND 16-DPSK MODULATED 
SIGNALS 


The introductory chapter presents an overview of the structure and performance 
of optimum receivers for M-ary signalling. The focus of this chapter is on DBDR’s. 
As mentioned in Chapter II, the symbol state generation based on the data state and 
the previous symbol state is the most important task performed by the encoder. 
Therefore in the receivers decoder, a similar algorithm must be performed in reversed 
Sense, that 1s, from two successive symbol states, the corresponding data state must be 
estimated. A generic form of the receiver in block diagram structure 1s shown in Figure 


Bal. 
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Figure 3.1 General receiver block diagram for DPSK signals. 


The algorithm for recovering the transmitted data is as follows 


1. Recover the symbol state from the received symbol using direct bit detection 
fearniques. 


2. Store the past symbol state oes 1, 2..., Vie 


3. Compare the present symbol state with the past symbol state to determine the 
incremental “ phase angle ~ between the two. 


* 


4. From Table 2 or 3 find the corresponding data state that gave rise 
to the “phase angle” between the two consecutive symbol states. 


For example, assuming the past symbol state is (.$,) 0001 and the present svmbol state 
is S$, (1001) from Figure 1.2.b, the angle between $, and S$, is 3a. From Table 3., the 
corresponding Gray Code representation for 3a is found to be (0010) , therefore the 
recovered data state will be D),. 

The decoder implementation can be depicted as a matrix structure in a manner 
similar to that shown for the transmitter encoder in Chapter II. For 8-DPSK and 
16-DPSK modulated signals, the decoder matrices are given in Figure 3.2 and Figure 
3.3 respectively. For example, assuming the same conditions as in the above 
paragraph, by using S$, and S, as a column and row entries for the matrix given in 


Figure 3.3, the recovered data state from the decoder is determined as D,, (0010). 


SYMBOL STATE 
oe eee ee 
SYMBOL 000 001 011 010 110 111 101 100 


s = 000 000 001 011 010 110 111 101 100 
<© _ 100 100 000 001 011 010 110 111 101 
<7? ~ 101 101 100 000 001 011 010 110 111 
<8 ~ 111 ~=«-111 101 100 000 001 011 010 110 
st ~ 101 101 111 101 100 000 001 011 010 
S2 = 010 O10 101 111 101 100 000 001 011 


5 = O11 O11 010 101 111 Pol 100 gooReut 
Es = 001 001 O11 O10R10L Ti Oreos oe 


Figure 3.2. Decoder matrix for 8-DPSK signals. 


B. COMPUTER GENERATION OF RECEIVER OUTPUTS 
Appendix C includes the results of the computer evaluation of the M-arv DPSK 
receiver. In addition to generating the decoder matrix in decimal form, the output for 


this program includes all possible combinations of the decoder outputs. 
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Figure 3.3. Decoder matrix for 16-DPSK 
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C. RECEIVER LOCAL OSCILLATOR PHASE ERRORS 

The introductory chapter summarizes the studies done on M-PSK DBDR 
structures, assuming that the receiver local oscillator has no phase error. In this 
research, a fixed receiver phase error has been assumed in the performance analyses 


conducted, and its effect on receiver BER is presented in Chapter IV and Chapter V. 
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IV. PERFORMANCE ANALYSES OF 8-PSK AND 8-DPSK DBDR’S WITH 
LOCAL OSCILLATOR PHASE ERROR 


In previous chapters, DBDR’s are introduced and their operational properties are 


explained. In this chapter, performance analyses for 8-PSK and 8-DPSK will be 
carried out in detail. 


A. 8-PSK DBDR PERFORMANCE ANALYSIS IN THE PRESENCE OF 
RECEIVER LOCAL OSCILLATOR PHASE ERROR 


The possible transmitted signals which are represented as vectors in Figure 1.2.a , 
can be mathematically expressed as 


s(t) = /2E/T, sin[ 2mfyt + 0.(t) le ie 3. 03 O0St st, (4.1) 
where 
ie 
Eee S-(tydt , i—ieens (4.2) 
) 
T, = symbol duration 
fy = carrier frequency (Hz.) 


Q(t) indicates the phase of the transmitted PSK signal, which for the 8-PSK case 
becomes 


a + (i-1) 7/4 i= 1,3,5,7 
0.(t) = O<StST (4.3) 
it/4-@ 1= 2,4,6,8 


The angle @ is variable as shown in Figure 1.2.a; however, the constraint 
0<Sa<7/4 must be observed. Equation 4.1 can be written by using a suitable 
orthonormal set of functions as 

s(= (0) VE, cos®,(t) + @,(t) VE, sin8,(t) (4.4) 


oT 


where @,(t) and @,(t) are the orthonormal function set which by applying the 
Gram-Schmitt orthonormalization procedure 
(Neto: pols 


can be shown to be 


©, (t) = A 2) sin 2nfyt 

p(t) = /2/T cos 2nf,t , Caan (4.5) 

The received signal depends upon the channel model and the transmitted signal 
as specified by Equation 4.1. In this research, an Additive White Gaussian Noise 
(AWGN) interference model has been assumed so that the received signal can be 


written as 
0) eS ean) (4.6) 


where n(t) 1s a sample function of a white Gaussian noise process having a constant 
(one sided) power spectral density function equal to Ny watts/Hz. The receiver 
processes r(t) over an interval of duration T, seconds and a decision is made at the end 
of the interval. This procedure continues in a repetitive fashion as time progresses. 
Assuming all signals are equally likely to be transmitted, the probability of a bit 


being correct for M-PSK signal transmission is written as 


M 
Pr{ bit correct.) = =e » Pr{ bit correct / s(t) transmitted } (4.7) 


1=1 


The calculations for the conditional probabilities (Pr{ bit correct / s(t) transmitted Ne 
i1=1,2,....M) have been determined by Myers and Bukofzer [Ref. 7: pp. 6-12], for the 
8-PSK case with no receiver local oscillator phase error. In the sequel, the same 
analysis methods are utilized to obtain the receiver BER in the presence of a local 
oscillator phase error. A local oscillator phase error is modeled by the angle € in the 
unmodulated receiver carriers @,(t) and @.(t) of Figure 1.6. These are denoted @, (t) 


and @,(t), namely 


38 


@,(t) = J2/T sin ( 2nft + €) 
P(t) = v2/T cos ( 2nft + €), OStST, (4.8) 


Assuming signal s.(t) was transmitted, this produces outputs Ty ; from the receiver's 
3 


correlators, where 
a 
tei = SUS(t) + a1 O, (t)dt , i=1,2..8, k=1,2 (4.9) 
0 


For example when s,(t) is transmitted, from Equation 4.4, Q.(t) becomes just @ and 


a VE, cos (a-€) + n, 


ry, = VE, sin (a-e) + n, (4.10) 
where 
oo 
n, =f n(t)@,(t)dt, k= 1,2 (4.11) 
0 


Observing that 


i 1 
E[n, ny] = E[ J n(t) @,(t) dt x J n(t) 9, (t) dt]= 0 (4.12) 
0 0 


This means that n, and n, are uncorrelated random variables. Since n, and n, are 
Gaussian random variables (having identical variances) [Ref. 7: p.7], they are also 
Statistically independent. 

The goal of the receiver of Figure 1.6 is to correctly recover the LB, the MB, and 


the RB, transmitted via the signal s.(t). By introducing the more compact notation 


Pr{ bit correct / s(t) } = Pr{ bit correct / s(t) transmitted } , 1=1,2..,8 (4.13) 
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the probabilities of such correct bit recovery are given in general by 


Pr{ tT ;> 0 / s.(t) Po 1=—12,8,7 
Pr{ LB correct / s(t) } = (4.14) 
PCG n seul 4-0 


Ri tector seer as En i= 1,4,5,8 
Pri MB correct / s(t) } = (4.15) 
Pr{ ty a =) 0, s(t) \ t= 2.3.64) 


Pr{ tape 107 S.(0) alii 12-324 
Pr{ RB correct / s(t) } = (4 16) 
mel tr, ,<0/ s(t) 1 = 5,6,1.8 


From Equations 4.8-4.16, we can write 


Pr{ LB correct / s,(t)} = Pr{ Tl 0 / s,(t) } 
= Pr{ n,> - VE, cos(a-e) / s,(t) } 


= 1-Q[ V2 E/N, cos(a-e) ] (4.17) 
where 
CO 
Q(x) = J [1/2] exp(-u7/2) du (4.18) 
x 
Similarly 


Pr{ RB correct / s,(t) } = Pr{ ae 0 / s,{t) } 
= Pr{ n,> - VE, sin(a-e) / s,(t) } 
1-Q{V2EJN, sin(a-e) } (4.19) 


In order to evaluate the conditional probability of the MB being correct, the following 


steps are required. From Equation 4.15 


Pr{ MB correct / s(t)} = Pr{r ee I, < >0/ s(t) } 
Pri(ry T> ;) ; (t) ;+ >) >0/ s(t) } (4.20) 
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By defining variables y, . and y, ; as 
: 7 


= 
| 


ee 5G 12; 


Yop = tart; (4.21) 


a gl 


then from Equation 4.9, when s,(t) is transmitted, 


i PREY [ cos (a-€) + sin (a-€) J+ n, +n, (4.22) 
i = 5 cos (@-€) - sin (a-€) J+ n, -n, (4.23) 


Clearly y, , and y, , are Gaussian random variables with statistical means given by 


iy) J E, [ cos (a-€) + sin (a-£) ] (4.24) 
Ey, ) = / E, [ cos (@-) - sin (a-£) ] (4.25) 


and variances given by 
Varf{ Yi = Var{ Yo = ING (4.26) 


By observing that 


eat ie E Ly) iD (Yo 1° E LY 1) pa Em (n, + n,) (nj- n,) } 
= E{(n,)*(n,)} = 0 (4.27) 


where n, and n, are defined in Equation 4.11. Since y, , and y,, are uncorrelated 
Gaussian random variables, they are also statistically independent. Therefore, 
Equation 4.20 becomes 


Pr{ MB correct / s,(t)} = Prf{ Yur You >0/ s(t) } 

Pr( Ly, ;>0/ s(t) ]and[y, ,>0/ s,(t)]}+ 

Pri [ yaa <0 / s,(t) ] and | oats O/s,(t)]}+ 

Prey) > O/ s(t) } Prey. >9/ s)(t) } + 

Eu one =0/ s)() 7 er Y24 07s (uy (4.28) 


(| 


+ 


I 


+ 
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From Equations 4.22 - 4.26, we have 


Pr{ y, > 0/ s(t)} = QLE{y yh / VNg] (4.29) 
Pr{ ¥5 > 0/7 s(t) = QO] Eloy W Nal (4.30) 


so that Equation 4.28 becomes 


Pr{ MB correct / s,(t)} = 1-Q{ ¥ E/N, [ cos(a-€) + sin(a-g] } - 
- Q{ V E/N, [ cos(a-€) - sin(a-€] } 
2 Oy AEM [ cos(@-&) - sin(a-s] } . 
Qt Vv E./Ng [ cos(a-€) + sin(a-€] } (4.31) 


Assuming now transmission of s.(t), a similar analysis procedure is performed 


but now Equations 4.14 - 4.16 become 


0,(t) = Th ae (4.32) 
tio = ¥ Pasi 2e) chen, 
Es ee cos(a+é) + n, : (4.33) 


where n, and n, are defined by Equation 4.11. From Equations 4.14 and 4.33 


Pr{ LB correct / s(t) } = Proxy 5-07 s(t) ) 
= Pr{n,> - VE, sin(at€) / s,(t) } 


1-Q{ J2EJN, sin(a+s) ] (4.34) 


and 


Pr{ RB correct / s,(t) } = Pri r,,>0/s.(t) } 


= Pr{n,> - VE, cos(a+8) / s,(t) } 
=1-Q[J2EJ/N, cos(ats) ] (4.35) 


For calculations involving the MB, one additional random variable needs to be defined, 


namely 


aR or, ae i= 1,2,..:.8 (4.36) 


yl 


Assuming transmission of s,(t), Equation 4.36 becomes 


ee (12.2 11,2 
=) E, eos (ate)tn, - J E. sin(at+e)-n, (4.37) 
where 
E {Y3 5} > Re cos (a+ &) - Es Sin(a@ - ©) =se {Yo 13 (4.38) 
Since 
Pr{ MB correct / s,(t)} = Pri (r, ,)° : (1, ,)°< 0 / s,(t) } (4.39) 


= Pr{ (Ty - ry >) (Ty + r, 5) <0 / s,(t) } 
by using Equations 4.21 and 4.36 
Pe B correct /s.(t)} = Prt Y1.2932< 0 / s.(t) } (4.40) 
Furthermore 


Ea (Y) 2° Ef Yy 23) (Y3 2° EA ¥323) aes (n, 1 n,) (n,- n.) } 
See) = 0 (4.41) 


so that y, , and y,., are statistically independent Gaussian random variables, and 


therefore Equation 4.40 becomes 


Bit Sib correct / s,(t)} = Prt ¥32 Yi2 < 0 / s.(t) } 

Pre [y,.79/ s)(t)] and[y;.>0/ sq(t)]}+ 

Pr{ [y¥,2<0/ s,(t)] and [y3 ,<0/ s,(t)]}+ 

= Prt yy 279 / so(t) } Pri y3 470 / so(t) + 

Pr{ y, 4<0/ s,(t)} Prt y3<0/ s,(t) } (4.42) 


+- 


+ 


By using a method similar to that used in the derivation of Equations 4.29 and 4.30, 


for y, 2 and ¥39> Equation 4.42 becomes 
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Piel B cormect 7 55(t) |’ uO VEG, [-cos(a@- 2)-  sint(@ae ae 
eh J ENG [ cos(a+ €) - sin(a+é] } 
= 2 Q(-y EMNG! cos(a+ 2) —> sin(@ee eine 
WO ERNE icos(@ ee) 5 cinema (4.43) 


Analvses can be carried out for probability of correct reception assuming every 
possible signal transmission. This is however not necessary because of inherent 
symmetries in the problem . In conditions where symmetries are not present, previous 
methodologies can be used to obtain results as follows. Assuming transmission of 


s,(t), Equations 4.3 and 4.9 become 


OF) it) 2s (4.44) 

3 = aes cos (@-E+ 7/2) + n, (4.45) 
= - ee sin (a@-£) + n, 

T,3 = ae sin (4-E+7/2) + n, (4.46) 


- JE, cos (a-£) + n, 
so that from Equation 4.14 


Pr{ LB correct / s(t) } = Pr{ r3< 0 / s,(t) } 
=- Eines ee sin (G-€) / s,(t) } 
= 1-Q[ v2 E,/Ng sin(@-) | (4.47) 


and from Equation 4.16 


Pr{ RB correct / s,(t) } = Pr{ T,3 < 0 / s,(t) } 
= Piaties a aie cos (-€) / s3(t) } 
= 1-Q) V2 E/N, cos(a+£) | (4.48) 


Assuming transmission of's,(t), Equations 4.3 and 4.4 become 
Q(t) = 1-@ (4.49) 


tp, = VE, cos (m-G-e) + n, (4.50) 
= -VE.cos(@+e) + n, 


| 


t,, = VE, sin (m-a-£) + n, (4.51) 
- JE, sin(a+e) + n, 


so that from Equation 4.14 


Pr{ LB correct / s,(t)} = Pr{ fe Os, {t) } 
—been, > 24/E. sin (a+ €)/ s,(t) } 
als Ory ues iNrmcosn(e: 16) | (4.52) 


and from Equation 4.16 
Piaweeeecorrect / s,(t)} = Pr{ T>4 > 0 / s,(t) } 


= Pr{ Wh Ee sin (a+ &) / S ,(t) } 
L- Ql V2 EJ/No sin (a+ €) | one 


Also from Equation 4.15 


Peehes correct /s,(t)} = Pr{ (T, 4) > ls ee UTES 8 (4.54) 
= Prt (ty g- Toy) Tg FPO / S40} 


Using Equations 4.21 and 4.36 vields 


(4.55) 
(4.56) 


ma tt 


mea 4 


Zio 


and using a method similar to that used in the derivation of Equation 4.27 results in 


EvvtBecorrect/s,(t)} = Prt Vie 4 >0/ s,(t) } 

Pr{ Ly, 4>0/ s(t) ]and[y, > 0 / sy(t)]}+ 

Pr{ [ ¥ g<0/ s4(t)] and [ y, 4<0/ s,(t)}}+ 

= Pr{ ata Se eet ¥2479 Hse) = 

Pi, Vie <0 / s4(t) } Pri Y24 <9 oe, Cire) 


alg 


+ 


which becomes 
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Pr{ MB correct / s,(t)} = 1-Q{ A PEIN, [ cos(a+€) - sin(a+6é) | } - 
-Q{ J2EJNp [ cos(a+e) + sin(a+e),J} + 
ae O) /2EJNg [ cos(a+eé)- sinfat+e)]}. 
.Q{ V2 E/N, [cos(ate) + sin(a+) ]} (4.58) 


For transmission of s,(t), the pertinent equations are 


8.(t) = m+ a (4.59) 

ts = y Eecosiitm 02s) sean (4.60) 
= - JE, cos (a-£) + n, 

To 5 = JE, sin (T+G-£) + n, (4.61) 


- Jie Gn (8-5) + 1, 
so that 


Pr{ MB correct / s.(t)} = Pr{ oa . (ere 0 / s(t) } (4.62) 
Pr{ (Ty 5 - > 5) (Ty 5 + Ty 5)? 9 | s(t) } 


which becomes 


Pr{ MB correct / s(t) } = Pr{ Vis Y25 >0/ s(t) } 

Prily, 57 0/ s(t) } and[y, 5>0/ s(t) ] 3+ 

Pr{ [ ie s(t) | and [ Y5<0/ S(t) p> 

Pr{ y, 5>0/ so(t) } Pr{ y,5>0/s,(t) }+ 

Pr{ ¥15 <0 [ sg(t) } Pri Y25 <9 / Sgt) } (4.57) 


+- 


+- 


by observing that y, , and y, are statistically independent random variables making 


it possible to apply the similar procedure shown above. The result becomes 


Pr{ MB correct / s.(t)} = 1-Q{ Ay eG [ cos(a-€) + sin(a-€) | } - 
-Q{ V2 E/N, [ cos(a-€) - sin(a-e) ]} + 
2 Oi aN, [ cos(a-£€) + sin(a-e)]}. 
-Q{ V2 E/Np [ cos(a-£) - sin(a-) J} (4.64) 
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As indicated before, because of the inherent symmetry in the problem, the 


following equations hold for conditional LB correct probabilities (see Figure 1.2.b ) 


Pri eB correct / s(t) } = Pr{ LB correct / s.(t) } (4.65) 
ae cemect / s5(t} } — Pr LB correct / S,(t) } (4.66) 
Pr{ LB correct / s,(t)} = Pr{ LB correct / s.(t) } (4.67) 
Pr{ LB correct / s,(t)} = Pr{ LB correct / S,(t) } (4.68) 


and similarlv for the RB 


Pr{ RB correct / s,(t) } = Pr{ RB correct / s.(t) } (4.69) 
Pr{ RB correct / s,(t) }} = Pr{ RB correct / s,(t) } (4.70) 
Pr{ RB correct / s,(t) } = Pr{ RB correct / s2(t) } (4.71) 
Pr{ RB correct / s,(t) } = Pr{ RB correct / s,(t) } (4.72) 


and finally for the MB 


Pr{ MB correct / s,(t)} = Pr{ MB correct / s,(t) } (4.75) 
Pr{ MB correct / s,(t) } = Pr{ MB correct / s,(t) } (4.76) 
Pr{ MB correct / s,(t) } = Pr{ MB correct / s,(t) } (4.77) 
Pr{ MB correct / s.(t) } = Pr{ MB correct / s.(t) } (4.78) 


By using now Equation 4.7 for the unconditional probability of the LB, MB and 
RB being correct, the above results yield 


Pr{ RB correct } = Pr{ LB correct } 


1- 1/4[ Qi V2 E/Np cos(a-€) J+ 

+ Q[ /2 E,/N, sin(a +s) ]+ 

+ Qi V2 E/N, cos(a+ £) ] 

+ QL V2 E/N, sin(a-€) J J (4.79) 


and 


Pr{ MB correct } = 1- 1/2[ Q{ V2 E,/Np [ cos(a-e) + sin(a@-€) ] } + 
+ Q{V2E/Npy [ cos(ate) + sin(at€) ]} + 
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+ Q{ VZEJNg [ cos(a-e) - sin(a-e) ]} + 
am ~/ ENING [ cos(a+) - sin(at+e) }} J+ 
oa J SEEN, [ cos(a@-&) + sin(a-£) ]}. 
OI </ OEINE [ cos(at+e) + sin(ate) ]} + 
+ Q (V2 E/N, [ cos(a-£) - sin(a-e) } } . 
no 2 E./Ng [ cos(a + €) - sin(at+s) }} ] (4.80) 


For 8-PSK,the Signal to Noise Ratio (SNR) has been defined by Equation 1.6, so 
that using Equations 4.79 - 4.80, the BER denoted PE can be expressed as a function 
of SNR, @ and € as follows 


PE= 1/3 [ Pr{ LB correct } + Pry MIB correct) seer, . os comect |} 
= 1-1/6 { Qi 2./SNR cos(a-e) ] + Q[ 2. /SNR sin(at+) J+ 

Qi 2 /SNR cos(a+e)] + Qf 2 /SNR sin(a-e) ] + 

+Q{2YSNR [cos(a-e) + sin(a-e) }}+ 

+ Q(2./ SNR [cos(at+e) + sin(ate)]}+ 

+Q{2/ SNR [cos(a-s) - sin(a-e) ]} + 

+Q{2/SNR [cos(at+6) - sin(ate)}}} + 

+ Q{2./ SNR [ cos(a-£) + sin(d-e) } }. 

. Q{2./ SNR [ cos(at+e) + sinf(ate)}]} + 

+ Q{2./ SNR [ cos(a-e) - sin(a-s) ] }. 

.Q{2Y SNR [ cos(até) - sin(ate) ]} (4.81) 


++ 


From Equation 4.81, the PE can be computed for various values of SNR, @, and 
€, with the results used to optimize the DBDR performance. For e= 0° (i. e. , no 
receiver local oscillator phase error), values of PE vs SNR have been plotted (see 
Figure 4.1) showing complete agreement with the results obtained by Myers and 
Bukofzer [Ref. 7: p. 11] Figure 4.1 includes also the PE for various values of £€, as a 
function of SNR. 


B. 8-DPSK DBDR PERFORMANCE ANALYSIS IN THE PRESENCE OF 
RECEIVER LOCAL OSCILLATOR PHASE ERROR 


The advantages and disadvantages of differential encoding and decoding have 
been previously explained in the introductory chapter. In this section, performance of 
the 8-DPSK DBDR 1s presented. As explained earlier the 8-DPSK decoder utilizes the 
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PHS. ERR.= 180 DEG 





=) 
—-4.0 1.0 7.0 iZco 18.0 LoL 
SNR (DB) 


Figure 4.1 PE Versus SNR of 8-PSK DBDR for Various Values of € (4= 2202): 
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LB, MB and RB, which are recovered via direct bit detection methods, and performs an 
inverse mapping from which the original data stream is delivered to the intended user 
Before proceeding with the performance calculations, the following notation will 


be introduced. The probabilities of a bit being correct or being in error will be denoted 


Ck = Pr{ k correct} , Ek = Pr{ k error } (4.82) 


respectively, where k= | for LB, k= 2 for MB, k= 3 for RB. By introducing the 
notation given by Equation 4.82, the analysis for the performance of the 8-DPSK 
DBDR can be expressed in more compact form. For example E2 means the 
probability that the MB is in error. 

In order to obtain the probability of a bit being correct (with differential 


decoding) we must first evaluate 


Pr{ RB correct } = 1/2 Pr{ RB correct / RB= 1} + 
+ 1/2 Pr{ RB correct / RB= 0} (4.83) 


where by further conditioning, we equivalently have 


loa 

Pr{ RB correct/RB=m} =) YY  Pr{ RB correct / RB=m,MB=i,LB=j }. 
i—0 si 

. Pr MB=1, LB— ae — Onl (4.84) 


Due to the assumption of equally likely probabilities of having 0’s and 1’s in the data, 
Pri MB=1, LB=j} = 1/4 1j= 0,1 (4.85) 


so that Equation 4.83 becomes 


ial 
Pr{ RBcorrect} = 1/8 ) YY Pr{ RB correct / RB=0,MB=i,LB=j } + 
i=0 j=0 
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l 
+ 1/8 ) YY Pr{ RB correct / RB=1,MB=i,LB=j } (4.86) 
ie 


The DBDRs recover from the received signal, each bit in the data state (for PSK 
modulation) or each bit in the symbol state (for DPSK modulation) independently. 
This property of the DBDR is an advantage for BER analyses. Because certain 
consecutive bit errors in the past and presently recovered symbol states do not 
necessarily result in errors in the recovered data states, these bit errors in successive 
svmbol states must be considered. That is, even though symbol errors may occur, the 
delivered data may stillbe correct. Therefore, these tvpe errors have to be accounted 
for as resulting in a correct reception and considered in the computation of the 
probability of a bit being correct. Figure 4.2 shows one example of the situation being 
discussed. Here the data state (0,0,0) is assumed to occur and the symbol state (0,0,0) 
is assumed to have been previously transmitted. This means that according to the 
encoder matrix shown in Chapter II, the svmbol state (0,0,0) must now be generated 
for transmission. Figure 4.2 further shows the receptions of the symbol states which 
even though, most contain errors, produce recovered data states in which the RB 1s 
correctly delivered. Observe that there are 32 such combinations which yield a correct 
QO for the RB. Therefore notation has to be introduced in order to account for error 
patterns of the above type, namely Pr{ RB correct, (k,lm) / RB=0,MB=1,LB=] }, 
where k,l,m= 0,1. The meaning of this expression is the probability of the RB 
correctly recovered given that the data state was (j,1,0) and the delivered data state is 
(k,l.m) (where m=0 in this case, otherwise the RB could not be correct due to the 
transmission of the RB=0). Therefore 


Pr{ RB correct / RB=0,MB=0,LB=0 } = 


ior 
= Y Pr{ RB correct,(k,l,0) / RB=0.MB=0,LB=0 } (4.87) 
k=0 1=0 


Since each element of the symbol state is independently recovered, we see from Figure 


4.2 that, for example 
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7HEEBBBAB BHAA AHAB AAA A 
CORRECT RECEPTION COMBINATIONS IN PRESENCE OF ERROR 
TERRE ABBB EAE AE AA AAA EBA BAB BREA R AAA AA HARARE FAA AAA ALHTTAAT TET 


INPUT PREVIOUS PRESENT PREVIOUS PRESENT OUTPUT 
TRANSMISSION TRANSMISSION RECEPTION RECEPTION 
0 0 0 0 20570 OC ae 0 


rrOOrFOORr ROOF Fr OORHPOORROORROORrRKROO 
Or OF OF OF OF OF OF OF OF OF OF OF OF OF OF OO 
OOOO RP RRR RrRrHOOOrF OOOOH HOR RRR ROOOO 
rrOOrrF OOF OF OOF OF FH OHPOOrF OR RROORROO 
Or OF OF OF OOF OR RP HOR RPOROOOrFRrOF OF ORO 
RRR RRR RRR RRR RR rH OOOOOOOOOOOO0O000 
rere rrr OOOO OOO OR PRR RRR rOOOOOO0OO 


rrr ODOOR RR RrHOOO OR RR RHOOOORRRRHOOO 
ODMOWOOODVDDWVODDDWVVDWVVDOVVOVOOVO OOOO CO0O0C00 


Figure 4.2. Example of Error Patterns for 8-DPSK System Where Data State, 
Previous Transmission, and Present Transmission are 000 with RB Correctly Recovered. 
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Pr{ RB correct,(0,0,0) / RB=0,MB=0,LB=0 } = Cl? C2? C3?+ 
+ Cl? C2* E3°+ Cl? E2? C3?+ Cl? E22 E32+ El? C2? C3?2+ 
+ El? C2? E3*+ El? E2? C3?+ El? E22 E32 (4.88) 


From Equations 4.79 and 4.82 - 4.86, it is clear that Cl= C3 and El= E3. 
Therefore, simplification of Equation 4.88 is possible, thus yielding 


Pr{ RB correct,(0,0,0) / RB=0,MB=0,LB=0 } = Cl1* C2?+ 2 Cl? C2? El?+ 
+ Cl* E2*+ 2Cl* El? E2?+ E1* C2?+ El* E2? (4.89) 


A similar expression must be obtained for other combinations of values of k and | as 
well as the condition that the RB is 1. This produces a total of 64 possibilities of 
received symbol states that result in a correct RB. A computer generated table was 
used ,in order to identify the equivalent error patterns for other cases involving MB= 1 
and LB=1. The appropriate listing are shown in Appendix C. Similar analyses must 
be carried out in order to determine the unconditional probabilities Pr{MB correct} as 


well as Pr{LB correct}. Such analyses demonstrate that due to the inherent symmetries 
Pr{ LB correct } = Pr{ MB correct } = Pr{ RB correct } = PE (4.90) 
so that using Equation 4.86 along with similar forms obtained for k= 1 and 1=1 yield 


PE = 1- 1/2[ 2C1* C2? + 9 Cl? C2? El*+ Cl? E2* + 3 C1* E2 C2+ 
+ 10 Cl? El? E2*+ 2 E14 C2? + El* E2? + 17 Cl? C2 El? E2 + 
+ C13 C22 El + 3 CIP El E2? + 2C1 El? C2? + 2C1 El? E2? + 
+ 3C2E1*E2 + 4C1 C2 El? E2 + 4Cl1° C2 El E2] (4.91) 


In Figures 4.4 and 4.5, PE (given by Equation 4.91 ) vs SNR have been plotted for 


various values of @ and € respectively. 


C. NUMERICAL RESULTS AND DISCUSSION OF RESULTS FOR 8-PSK 
AND 8-DPSK DBDR’S WITH RECEIVER LOCAL OSCILLATOR PHASE 
ERROR. 


Previous sections of this chapter have focused on the performance of 8-PSK and 
8-DPSK DBDRs operating in the presence of a receiver local oscillator phase error. 


Without any phase error, we have shown the performance plots of 8-PSK and 8-DPSK 
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Figure 4.3 PE vs SNR Comparison Plot for 8-PSK and 8-DPSK DBDR 
(= 22.5° , €= 0°). 
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Figure 4.4 PE Versus SNR of 8-DPSK DBDR for Various Values of ¢ (e= 0°), 
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Figure 4.5 PE Versus SNR of 8-DPSK DBDR for Various Values of € (@= 22.5°). 
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DBDRs on Figure 4.3, thus demonstrating the performance penalty of the latter 
receiver. While usage of 8-DPSK provides good bandwidth efficiency (same as 8-PSK) 
and relatively good noise immunity, the BER of 8-PSK is lower than that of 8-DPSK 
for €= 0°, as expected 

The value of using 8-DPSK versus 8-PSK can be appreciated more fully when 
comparing the performance plots of Figures 4.1 and 4.5. These show that if the 
receiver has a local oscillator phase error restricted to 0S€566.5°, for similar SNR 
values, the BER of the 8-DPSK receiver is lower than that of the 8-PSK receiver. 
Since such phase errors can be tvpically expected to be small, this is a significant 
advantage. However it must be pointed out that for phase errors of more than 10°, 
both receivers exhibit severely degraded performance that would make both svstems 
unsuitable for operation in a practical svstem. The 8-DPSK receiver does not exhibit 
complete insensitivity to phase errors as encountered in phase measurement receivers. 
Only when a complete phase revesal takes place (1.e., €= 180°) is the performance of 
the 8-DPSK DBDR unaffected, while the 8-PSK receiver produces a nearly useless 
output. 

The angle @ set to 22.5° is optimum in minimum BER sense for the 8-PSK 
DBDR. The variations in PE versus @ are shown in Figure 4.6, for ¢= 0° for 8-PSK, 
and a similar plot of PE versus @ is plotted for various values of SNR and €= 0° on 
Figure 4.7 for the 8-DSPK case. 


a 





15.0 18.0 aieo 24.0 Lio 30.0 
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Figure 4.6 PE Versus a of 8-PSK DBDR for Various Values of SNR (c= 0°). 
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Figure 4.7 PE Versus « of 8-DPSK DBDR for Various Values of SNR (€= 0°). 
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V. PERFORMANCE ANALYSES OF 16-PSK AND 16-DPSK DBDR’S 
WITH LOCAL OSCILLATOR PHASE ERROR 


In this chapter, DBDRs will be analyzed for 16-PSK and 16-DPSK 
communication systems operating in the presence of a receiver local oscillator phase 
error. The analysis procedure will be similar to that carried out in the previous chapter 
for 8-PSK and 8-DPSK systems. 

The performance of 16-PSK DBDR’s has been investigated and the methodology 
for carrying out the performance analyses have been well-documented by Bukofzer 
[Ref. 8], In this research, these analyses will be extended to the 16-DPSK system as 
well as the 16-PSK system modified to operate in the presence of a receiver !ocal 


oscillator phase error. 


A. 16-PSK DBDR PERFORMANCE ANALYSIS IN THE PRESENCE OF 
RECEIVER LOCAL OSCILLATOR PHASE ERROR 


The possible transmitted signals, which are represented as vectors in Figure 1.2.b, 
are mathematically given by Equation 4.1 with 1= 1,2,3,...,16. The same notation 
introduced in Chapter IV will be used to represent the signals using the orthonormal 


set given by Equation 4.4. Therefore the transmitted signals are expressed as 
s(t)= @,(t) VE, cos6.(t) + ,(t) JE, sin@,(t) (5.2) 
where 1= 1,2,...,16, OS tao and 
a+ (i-1) 7/8 i= 1,5,9,13 
B + (1-2) 7/8 i= 2,6,10,14 
0.(t) = 4 (it 1) 7/8 -B 1=3,7,11,15 (5.3) 


17/8 -@ i= 4,8,12,16 


with @ and B as shown on Figure 1.2.b 
Assuming all signals are equally likely to be transmitted, for the 16-PSK case 


Equation 4.7 can be written as 
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16 
Pr{ bit correct }} = 1/16} Pr{ bit correct / s(t) transmitted } (5.4) 
1= 1 


The notation to denote the four bits that make up a data state has been introduced in 


Chapter I. Also the notation introduced by Equation 4.13 is used here, namely 
Pr{ bit correct / s(t) } = Pr{ bit correct / s(t) transmitted } i= 1,2..,16 (Ses), 
From Figure 1.2.b, the conditional probabilities for the correct reception of the LB, 


MBL, MBR, RB are given by 


Pr{ LB correct / s.(t) } (5.6) 


Pr{ tT) > 0/ s(t) } i= 1,2,3,4,13,14,15,16 
Pr{ ty jes s(t) } i= 5,6,7,8,9,10,11,12 


Se ON S(t) = 12,7; Spo momo 16 


Pr{ MBL correct / s(t) } (Sa) 


Pr{ ry jr seema0ues.( ts) i= 3,4,5,6,11,12,13,14 


Pr{ MBR correct / s(t) } = 


Pr{ Pie / s(t) } i= 1,4,5,8,9,12,13,16 
(5.8) 


Prt V7}, ~0'/ s(t) } 1= 2,3,6,7,10,11,14,15 
Pr{ r,.>0/ s(t) } i= 1,2,3,4,5,6,7,8 
Pr{ RB correct / s(t)} = (79) 
Pr{ r,,<0/ s(t) } 1= 9,10,11,12,13,14,15,16 
where the r, , are given by Equation 4.9 and 
Von = (0 Pty? 2-4 (pty?) 15 1,2,3,-16 (5.10) 
The analyses methods which have been given in Chapter IV, can be applied to 


determining the conditional error probabilities for LB, MBL, RB. For MBR, the 


analysis procedure which has been used by Bukofzer [Ref. 8], will be utilized. In order 
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to define the conditional probability sets of Equation 5.8, additional steps must be 
completed as demonstrated in ihe sequel. If s(t) transmitted, the MBR 1s correctly 
recovered if for i= 1,4,5,8,9,12,13,16, V7, 20 andif for t= 2,3:G;7e00;) tea ietey ienes 
Even though the r nj are conditional Gaussian random variables, because of the highly 
nonlinear form of Equation 5.10, a more productive approach to the determination of 
the conditional p.d.f of V7, has been introduced (see Reference 8 for details) rather 
than using transformations of Gaussian random variables. By introducing the 


following notation, namely 


r) ;= component of r(t) along @,(t) = Ecos n (Sas 
r,.= component of r(t) along @,(t) = sin n (5942) 
it can be demonstrated that 
ry = C* cos 4n (5.13) 
where 
€ = length of the vector = V(r, 3° + (r5 ,)° (5.14) 
nN = vector angle that with respect to the positive x axis (5:15) 


= arctan (Tr, / Ty 3) 


Since €*> 0, it is clear that the test involving Vy}, being greater than or less than Zero 
is noW equivalent to the test involving cos4n in the same manner. This implies that 
the p.d.f. of mn must be obtained which, as shown in Reference 8, leads to integrals that 
cannot be evaluated in closed form and therefore must be evaluated via numerical 
integration. 

By utilizing the analysis procedure from Reference 8 and a methodology 
introduced in the previous chapter for the 8-PSK case, the conditional probabilities for 
the correct reception of a bit in the presence of receiver local oscillator phase error is 


given in shorthand form by the following expressions 


Pr{ LB correct / s(t)} = 1- Q{ /2 E/N, f(a,B,e)}, i=1,2,3...16 (5.16) 


where 
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cos (G-€) Tee 


cos (B-€) i= 2,10 
sin (B+ 5) i= 3,11 
f(a.B,e) = ¢ sin (a+e) i=4,12 (5.17) 
sin (B-e) i=5,13 
sin (@-£) 1=6,14 
cos (B+ ¢e) i=7,15 
cos (a+ €) 1= 8,16 


and @ and $B are shown on Figure 1.2.b. Furthermore 


Pr{ MBL correct / s(t) } = Q{-v E./N, (cos y+sin y) . Q{-VE,/Ny (cos y-sin y) } + 
Oy Ee Nancos y-siny) ). Ol EJ Na(cos y+ sin ¥)} (5.18) 


where 


y= on i=1,5,9,13 (5.19) 
(B-e) i=2,6,10,14 


Or 


Piqeieie correct / s(t) } = 1 - Q{-./E/N, (cos y+sin y) . Q{-VE./N, (cos y-sin y) } + 
+ Q{(VE,/Ny (cos y+sin y)} . Q{VE/Ng (cos y-sin y)} (5.20) 


where 


es eee (5.21) 
(a-€) i1=4,8,12,16 


For the calculations involving the conditional probability for the MBR, the joint 


p.d.f. of the random variables n and £ has been determined by Bukofzer [Ref. 8] 
namely 


ens, (er) = iL} fF tare (L cosH , L sinH) + (622) 


+ |LIf ee iecoshl = basi atil 
pin aus: 
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which must be integrated to produce the marginal p.d.f. of n. This results in 
IVs, (/s.) = 12 ew + fit J dE JEN, - (3:23) 
. { exp[ E/N) cos” (H-§,)] cos(H-&) . Ql -/2 E/Ny (cos(H-§)] } 


where 
6. = O(t)-¢ On (9.24) 


where 0.(t) is given by Equation 5.3 Therefore, the probability that V7,>0 involves 
computing the probability that cos4n>0. Observe that cos4n>0 for -7/8<n<17/8, 
37/8 <n < 50/8, 72/8< Nn < 90/8, and 11t/8<n<13n/8. Therefore the calculation of 
the conditional probability of the MBR being correct given that s,(t) was transmitted 
can be calculated by integrating the conditional p.d.f of n over the appropriate 


regions, namely 


Pr{ MBR correct / Sry) Pr{ cos4n>0 / oO = 


1/8 50/8 97/8 
= J fyjs, Aisi] SH + J fy. (His) dH + f fy), (H/s,) dH + 
-10/8 37/8 77/8 
13/8 
ae n/s, [H/s,] dH 
1172/8 
7/8 
=Pb i ae n/s, [H /s] + £ ns, (H+ 1/2) / s)] + fn/s, (H+) /s,] + (S253 
-7/8 


+ fyjs, (H+ 32/2) / s,]} dH 


By utilizing the appropriate form for the conditional p.d.f. of n when s(t) is assumed 


transmitted, results in 


Pr{ MBR correct / s(t) } = 1/2 exp(-E./Ny) + 
7/8 
+ J E/tNy J { exp[- E./No sin?(H-2)] cos(H-A) . {QI-./2 E./Np (cos(H-)) - 
-1'8 
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- Qi/2 EAN  (cos(i+h)] } + { exp[-E./No cos(H-2)] sin(H-d) Q[-/2 E/N, (sin(H-A)] - 
- QV2 EIN, (sin(H +4) ]} } dH (5.26) 


where 


ate  i=4,8,12,16 

h =la-e i= 1,5,9,13 (5.27) 
B-e i=2,6,10,14 
B+e  i=3,7,11,15 


Finally, it can be seen that 
Pr{ RB correct / s(t) } = Pr{ LB correct / s,(t) } (5.28) 


so that by using Equations 5.4 and 5.16-5.28 for correct reception probabilities of the 
individual bits, the overall BER can be obtained by using 


Pa, Oi@in error} = | - Pr{ bit correct } (5:29) 
This results in the (unconditional) probabilities 


feiyeteoein error } = Pr{ RB in error } 
= 1/8 { Q V2 E/N, cos(a-e)] + Ql /2 E/N, cos(B-e)] + 
Oly 2 EY Nec) el Oly 2 EN. sina s)] + 
+ Qi V2 E/N, sin(B-s)] + QLV2 E/N, cos(B+8)]+ 
PeOLy QE Ne cos eae Oly 2 ENE sing + s)) (5.30) 


and 


Pr{ MBL in error} = 1/2 + 1/4[ Q {-VE,/N, [cos(a-e) + sin(a-e)]} + 
+ Q {- VEIN, [cos(a+£) - sin(at+ €)]} + Q {-VE,/Nj [cos(B-2) + sin(B-e)]} + 
+ Q (-VE,/No [cos(B-e) - sin(B-e)]} - Q {-VE,/Ng [cos(B+ €) + sin(B + €))} - 
- Q {-VE,/N, [cos(B + €) - sin(B + €)]} - Q {-VEJNg [cos(at+e) + sin(a+)]} - 
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e {-VEJ/N, [cos(a +€) - sin(at+ €)}} ] + 1/2[ Q {-VE,/Ny [cos(a+e) + sin(a+ €)]} . 
YEN, [cos(a-€) - sin(a-£)]} + Q {-VE/Ng [cos(B+e) + sin(B+e)}} . 
'E /No [cos(B + €) - sin(B + €)]} - Q — [cos(a-€) + sin(a-s)]}} . 
50) 7 ER [Nog [cos(a-€) - sin(a@-€)]} - Q {-VE./Np [cos(B-e) +  sin(p-e)]} . 
. Q {-VE./Ng [cos(B-€) - sin(B-€)]} ] (5.31) 


For the MBR the result appears in integral form. That is 


7/8 
Pr( MBR in error} = 1-1/2. E/2N, § { exp{-E,/N, sin*[H-(a-e)]} cos[H-(a-e)]. 
-7/8 
{ 1/2 - Ql-VEJ/N, cos(H-(a-e))]} + exp{-E,/N, cos*!H./a-e)]} sin[H-(a-e)] 
.{ 1/2 - Q-VEJ/N, sin(H-(a-e))]} + exp{-E,/N, sin*[H-(B-e)]} cos[H-(B-e)] 
( 1/2 - Q- E/N, cos(H-(B-e))]} + exp{-E,/N, cos*[H-(B-)]} sinfH-(B-e)] 

. {1/2 - Q-VEJNg sin(H-(B-€))]} - exp{-E,/N, sin*[H +(B + £)]}} cos[H + (B+ )] 

.{ 1/22 0[2 JEN, cos(H + ($+ €))}} - “ak ,/No cos*{fH +(BP+e)}} sinfH +(B +e) 
{ 1/2 - Q[- nee Ne sin(H +(a+e))]} + exp{- E No sin*[H + (a+¢)]} cos[|H+(a+6)] 
{ 1/2 - QI- ae Ne cos(H+(a+6)]} + exp{-E./N, cos 7H +(a+¢)]} sinfH+(a+ 6) 
( 1/2 - Ql-VE./N, sin(H+(a+))]} } dH (5.32) 


The overall BER can now be determined from Equations 5.30,5.31 and 5.32 by 


evaluating the final result for PE, the bit error probability, namely 


PE = 1/4{ Pr{ LB in error} + Pri MBL imeuneme, 
Pr, MBR in error } ebm RB im crer se, (5.33) 


It can be seen that Equation 5.33 consists of many terms that by themselves do not 
yield as much insight as the computed values of PE presented in plotted form. 
Therefore plots of BER vs. SNR for various values of @ and &, are given as Figures 5.1 
and 5.2. While a@ and B can be set independently, from an optimization point of view, 
a+ P=7/4 has been shown to produce the lowest BER for given values of SNR when 
€= (0° [Ref. 8] 
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Figure 5.1 PE Versus SNR for 16-PSK DBDR (a=11.25°, ¢=0°). 


67 


10° 


=| 


10 


Figure 5.2. PE Versus SNR for 16-PSK DBDR for Various Values of ¢ (a= 11.25%). 
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B. 16-DPSK DBDR PERFORMANCE ANALYSIS IN THE PRESENCE OF 
RECEIVER LOCAL OSCILLATOR PHASE ERROR 


Similar analyses on DBDR performance have been carried out for 16-DPSK (see 


Chapter I[V.B), so that the notation introduced previously will be used here. That is, 
ees COrrect } 5 EK = Pr { k in error} (5.34) 
where now k=1 for LB, k= 2 for MBL, k=3 for MBR and k= 4 for RB. 
In order to evaluate Pr { LB correct } with differential decoding, we have to 


evaluate 


Pr, LB correct }= 1/2 [ Pr{ LB correct / LB=0} + 


+ { Pr{ LB correct / LB=1 } ] (3.35) 
where 
Pr{ LB correct / LB=0}= (5.36) 
ea ] 
=) ¥ YD Pr{ LB correct / LB=0, MBL=i, MBR=j, RB=k} 
1=0 j=0 k=0 


x Prf MBL=i, MBR=j, RB=k} 


where for i,j, k= 0, 1 


Prf MBL=i, MBR=j, RB=k}= 1/8 (5.37) 


In Chapter [V, it has been explained that with DPSK modulation in conjuction 
with bit direct detection methods, certain bit errors in the past and presently recovered 
symbol states do not result in errors in the recovered data states. Such bit errors in 
successive symbol states must be considered, because, even though symbol errors may 
Seemmetie Gelmered data may still be correct. {herefore, these errors have to be 
accounted for and considered in the computation of the probability of a bit being 
correct. Figure 5.3 shows one example of the situation being discussed. That is, the 


data state (0,0,0,0) is assumed to occur and the symbol state (0,0,0,0) is assumed to 
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TBR EER RRA Brie AERA AEE EEE ES AE TESA EEE EEE SEES PSE TP TIE TPE cae ae ECT ESE 


CORRECT RECEPTION COMBINATIONS IN PRESENCE OF ERROR 


(HARARE AAA BARBARA AAA BHA AAA AAA HAA BE Ara HAAR ETETE 


OUTPUT 


PRESENT 


RECEPTION 


RECEPTION 


PRESENT PREVIOUS 


PREVIOUS 
TRANSMISSION TRANSMISSION 


TNeU 


: DCO OVOVOSVSSTOOVODVSOOSCODODODOODODOOODOODCOCOD ODDO OD DOOD DOOD COCO COCO OOK O00K0 
ge ee Se? OS OOO NAAN AN A ANNA AA AAA ANAS OS ODOOCOOOOO OCS OC Nd sate Ia asl 
‘ CTO COO SO DODO DOODDODOODVCOCO COCO OC CO OOO Ott Ht ttt HAA Att det ddd tetedeted 
’ OCOOSSSSOSSTOSVVDVVOVVCODVOODDVODODDODVODDCOO DDO OOD DODO O OOD ODO OCOD OOOO COCOOO0D0 


AFOADADADADAD AOD ADDO ADDO O HAO AAA AAAOAODAODOOOHOHOHOHADDOOOODOOOHHAAeoaae 
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ODVDVDVOV OC FHA HA AHH AAA IAD ADD OAD AOQAO A AAAAO AAAHOOAHOOOOOHOHOHOHOHOHOHO 


Geos One 
0 


FAOADAOAOAOAO AOA OAHOAO ADAH OAOTHO AO HOHOHOHOHOHOHOHOHOHOHOHONHOHOHOH 
OAADD AAD O AAD DAADOAADO AADCOAHOOAAOO AHO OHAAOQOOAHHOOAHOOHAHOOHHOOHHOO 
COOSAAAADOD DOA AAADOO ONAN AAD OOO HAAAAD OOO HAAAHOOOO NH HAAHOOOO HAHA AOOOO-e 
COSSS SDAA AAAAAAD OOOO OOO NM Ht AH AAMAAHOODOD OOOOH HAHAHA AAO O OO OOOO FAHetoded 


0000 


0/0500 


0000 


0000 
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Previous Transmission, 


Figure 3 


and Present Transnussion are 0000 with RB Correctiv Recovered. 
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have been previously transmitted. This means that according to the encoder matrix 
shown in Chapter II, the symbol state (0,0,0,0) must now be generated for 
transmission. Figure 5.3 shows the receptions of the symbol states which produce 
recovered data states in which the RB is correctly delivered. Observe that there are 128 
such combinations which yield 0 for the RB. Therefore notation has to be introduced 
to account for error patterns of the above type, namely Pr{ RB correct, (g,h,p,q) / 
LB=m,MBL=1,MBR=],RB=k 3}, where g,h,p,q= 0,1. The meaning of this 
expression is the probability of the RB being correctly recovered when the data state is 
(m,j,1,k) and the delivered data state is (g,h,p,q) where q=0 in this case (otherwise the 


RB could not be correct due to the transmission of the RB=0). Therefore 


Pr{ RB correct / LB=0, MBL=0 “3R=0, RB=0} = 
Ly ie 


= YY Pr{ RB correct, (g,h,p,0) / LB=0, MBL=0, MBR=0, RB=0 } (5.38) 


s=0 h=0 p=0 


Since each element of the symbol state is independently recovered, we see from Figure 


5.3 that 


Pr{ RB correct,(0,0,0,0) / LB=0,MBL=0,MBR=0,RB=0 } = Cl? C2? C3? C4*+ 
+ Cl? C2* B3* C4? Cl* G2? Bs? E42 G12? @aeieqs Classe serene 


+ Cl? B2? £3? Bd" El? + €2* C3? C47 Bi C2@s bd Bene. ba ere 


+ El? C2? £3? Bd? Bi? B2* C3? C441 2” Eads Se 
+ Cl? C2? C3? E4? (5.39) 


rlowever, from Equations 5.34, and 5.30, it is shown that Cl= C4 and El= E4. 
Therefore: 


Pr{ RB correct, (0,0,0,0) / LB=0, MBL=0, MBR=0, RB=0 }= 

Cl? C22 C3*+ 2Cl C2* G37aEl* + C2 tee” bie eee 

2 C14 C2? E32 + 2Cl* C2? El* E3*+ Cl* C32 E22 + El? E22 E3? 

2 Cl? C3? El? E2? + C2? El? E2? E32 + Cl? El? E2? E37+ 

C3? E1* E2* + Ci? C3? E14 (5.40) 


tf + I 


-{- 


ee 


‘A similar expression must be obtained for other combinations of values of g,h and p as 


well as the condition that RB is 1. This produces a total of 256 possibilities of received 


svmbol states that result in a correct RB. Computer generated outputs were used to 


identify the equivalent error patterns for other cases involving LB=1, MBL=1 and 


MBR=1. These results are shown in Appendix C. Similar analyses must be carried 


out in order to determine the unconditional probabilities Pr{LB correct}, Pr{MBL 


correct} 


as well as Pr{MBR correct}. Such analyses demonstrate due to inherent 


symmetries that 


Pr{LB correct} = Pr{MBL correct} = Pr{ MBR correct} = Pr{RB correct} = PE (5.41) 


so that using Equation 5.38 along with similar forms obtained for k=» i and 1=1 yield 


PE = 


fi 
z 
+ 
+ 
+ 
+ 
+ 
i 
2 
is 
+ 
+ 
+ 


ete ©) Gi el” C22 G3. E3 £93°Cl* C2? C3 E2+ 

2C1* C2 C3 E2 E3 + 6 Cl? C2* C32 El + 4 Cl? C2? C3 El E3+ 

Piel G @s- Gs El E2 + 8 Gl G2 @3 G3 El E2 + 6 Cl’ C2* El E3°+ 
6 C13 C2? C3* El + Cl* C3 E2? E3 + 4 C1? C3 El E2* E3+ 

peri Es- + Cl* C2> E3* + 3 Gl 2 E?2 E3* + 2 C2? C3* El*+ 

@l eo Ee + 10 C1? C22 C3 El* + 10 Cl? C2? El* E3* + 

10 Cl* C3? El* E2* + 8 Cl? El? E2* E3* + 18 Cl* C2 C3? El? E27+ 

18 Cl* C2 El? E2 E3* + 6 Cl? El E2? E3? + 12 Cl? C2 El E2 E3*+ 

12 Cl? C2 C3 El? E2 E3 + 6 Cl? C2* C3 El? E3 + 5 Cl? C3 El? E2? E3+ 
6 Cl C2? C32 El? + 6 Cl C3* El? E2? + 6 Cl C2? El? E27+ C2? El? E32 
emne? @s- le Er? Ger El B2Es- agG) El” E2, E3-+ 
me? Geel E2 Es Cl Cs El ee?, E3 4 Cl C2* C3 El’ E3+ 
C32 El* E2?+ 2 El* E2? E3* + C2? C3 El* E3 + C3 El* E2* E3+ 

3 C2 El* E2 E3* + 2 C2 C3 El* E2 E3 + 2 C2 C3? El* E2} (5.42) 


In Figure 5.4, PE (given by Equation 5.42) is plotted for various values of @. Also 


Figure 5.5 includes BER for various values of € when @= 11.25° 
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Figure 5.4 PE Versus SNR (dB) for 16-DPSK for Various Values of a (€= 0°). 
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Figure 5.5 PE Versus SNR (dB) for 16-DPSK for Various Values of ¢ (a= 11.25°). 
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C. NUMERICAL RESULTS AND DISCUSSION OF RESULTS FOR 16-PSK 
AND 16-DPSK DBDR’S WITH RECEIVER LOCAL OSCILLATOR PHASE 
ERROR 


In this chapter 16-PSK and 16-DPSK DBDRs have been analyzed when 
Operating in the presence of a receiver local oscillator phase error. The performance 
results which have been obtained for the 16-ary case, are similar to those obtain for 
the 8-ary case. Hence, instead of repeating the conclusions similar to these presented 
for the 8-ary case, conclusions that are unique to the 16-ary case will be presented in 
this part. Main differences between 8-ary and 16-ary cases which are observed via 
analyses and numerical evaluations are 


1. The optimum angles in minimum BER sense, @ and B between the symbols 
becomes 11.25° and 33.75° respectively for 16- PSK — case, as well as for 
16-DSPK case, with a+ B= 45° 


2. The usage of 16-PSK provides better performance than 16-DPSK for € S$ 33.5°, 
however the actual BER values are too high for use in a practical system when 
€ greater than a few degrees. 


3. The 16-DPSK system performs better than the 16-PSK system for €> 33.5°. 
However again the actual values of PE under these circumstances is too high 
for a practical system. 


4. The 16-DPSK system is insensitive to phase errors of 7 radians. 
Performance comparisons between the 16-PSK and 16-DPSK DBDRs are plotted 
on Figure 5.6, and BERs are plotted on Figures 5.7 and 5.8 respectively for various 
values of @ with e= 0° 
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Figure 5.6 PE Versus SNR Comparison Plot for 16-PSK and 16-DPSK DBDR 
(a= 11.25°, €= 0°). 
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Figure 5.7. PE Versus @ for 16-PSK DBDR for Various Values of SNR (€= 0°). 
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Figure 5.8 PE Versus a for 16-DPSK DBDR for Various Values of SNR (c= 0°). 
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VI. CONCLUSIONS 


In this research the use of DBDRs had been analyzed for M-PSK and MDPSK 
communication systems. The main assumptions in the analyses are a Gray coding 
mapping for the symbols and an Additive White Gaussian Noise channel model. It is 
further assumed that the receiver has a known phase error in its local oscillator and its 
effect on receiver performance is analyzed for the the special cases in which M=8 and 
M = 16. 

In addition to the analytical results, graphical results are presented in the form of 
plots of variations in PE for various values of phase angle a, SNR, and local oscillator 
phase error €. 

Optimum conventional M-PSK and M-DPSK receivers are designed to detect 
symbols while DBDRs are designed to detect bits. However, local oscillator phase 
errors degrade the performance of M-PSK and M-DPSK conventional receivers. Part 
of the work undertaken was intended to investigate the effect of such phase errors on 
DBDR’s. The variations of PE versus € are plotted in Figure 6.1 for 8-PSK and 
8-DPSK, and in Figure 6.2 for 16-PSK and 16-DPSK. If € (receiver local oscillator 
phase error) can be constrained to a small interval, namely 0° S¢< 30°, the M-PSK 
system provides better performance than the M-DPSK system. For values of € outside 
this range, both systems have high error rates that make the receivers unsuitable for 
practical applications. However, for complete phase reversal, i.e., €= 180°, the M- 
DPSK system suffers no performance penalty whereas the M-PSK system can not 
properly operate. 

The performance degradation due to a phase error has been analyzed also from 
a different point of view. That is, when € # 0°, there is an increase in PE for both the 


M-PSK and M-DPSK receivers. The ratio of the phase error is given by 


PER = PE, | PE, (6.1) 
where 

PE, = Probability of bit error @ € = 0° 

PE, = Probability of bit error @ € = 0° 
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PER has been analyzed under various conditions for M-PSK and M-DPSK. The 
results have been plotted in Figures 6.3and 6.5 for M=8, and in Figures 6.6-6.8 for 
M=16, as a function of SNR for different values of &. The plots demonstrate 
significant differences in the performance degradation of the M=8 versus M = 16 case. 
However, for high SNR values, these differences tend to vanish. 

The analysis results lead to the following conclusions 


1. By using DBDRs, the error probability of each bit in the symbol state can be 
calculated independently and more importantly these error probabilities can be 
varied by changing the angle between symbols. 


2. For optimum BER, @ should be set to 22.5° for 8-PSK and 8-DPSK systems, 
while 16-PSK and 16-DPSK systems require @= 11.25° and B+a=45°. 


The solution for the phase ambiguity problem cannot be achieved by using M- 
DPSK DBDR’s, unless the phase ambiguity is 180°. Therefore, if the perfect 
coherency conditions cannot be achieved, the M-PSK system should be used 
Whenever DBDR’s are used. 


G2 


4. When e= 0°, the DPSK system exhibits good performance. However, if the 
€= 0° condition can be guaranteed there is no need to use the differential 
encoding scheme. 
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Figure 6.1 PE Versus € for 8-PSK and 8-DPSK DBDR (SNR= 15 dB and @= 22.5%). 
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Figure 6.2 PE Versus € for 16-PSK and 16-DPSK DBDR (SNR= _ 15 dB and @= 11.25%). 
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Figure 6.3 PER (e€=0°/e=5°) Versus SNR for 8-PSK 
and 8-DPSK DBDR’s (a= 22.5°). 
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Figure 6.4 PER (e€=0°'e=22.5°) Versus SNR for 8-PSK 
and 8-DPSK DBDR’s (a@= 22.5%). 
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Figure 6.6 PER (€=0°'e= 5°) Versus SNR for 16-PSK 
and 16-DPSK DBDR’s (@= 11.25°). 
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APPENDIX A 
DBDR INPUT/OUTPUT GENERATION PROGRAM (BELIZ) 


1. USER GUIDE 
This program generates all the possible input / output combinations in an M- 
DPSK communication system. The program runs interactively and sends outputs to 
the files 
BOLE FLOsE00! 
FILE RCDATA 
3. FILE XMDATA 
A) allem 
peewee ERR. 


These files are used for information transfer between subroutines but are not required 


to 


for program execution. It is possible to generate input / output combinations for 
8-DPSK and 16-DPSK by using this program, while for 32-DPSK or 64-DPSK 
svstems, dimensions and format statements must be changed accordingly, and a larger 
character variable set in the subroutine XMTR must be created. The program accepts 
only the integer equivalents of vour data or svmbol states. An example run of the 


program 1s shown below 


run beliz 
MABCULION BEGINS... 


PEBRSe SELECT ONE OF THEM 
3 he YOU RAVE TRANSMITIER 
2) fe YOU HAVE RECEIVER 


A 
PereotsSeEECT ONE OF THEM 
+3} Pec OUneorolEW S-DPSK 
2) ly OUR SYSTEM Lo-DPSK 


1 
ENTER THE INTEGER EQUIVALENTS OF SYMBOL STATES 
ithe ORDER OF COUNTER CLOCK WISE 

ENTER THE STATE # 1 


Elm@ek THE STATE & 2 
ENTER THE STATE # 3 
ENTER THE STATE # 4 
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ENTER THE STATE # 5 


2 

ENTER THE STATE # 6 
0 

ENTER THE STATE # 7 
4, 


ENTER THE STATE # 8 


6 
ENTER THE INTEGER EQUIVALENTS OF THE 


PHASE ANGLES OF YOUR CODING S7SrEM 
1*ALPHA= ? 


0 

2*ALPHA= ? 
1 

3*ALPHA= ? 
3 

4X*ALPHA= ? 
2 

5*ALPHA= ? 
6 

6*APLHA= ? 
7 

7*ALPHA= ? 
5 

8*ALPHA= ? 


4 
PLEASE ENTER THE SELECTION CORRECTLY 
ENTER 1 IF YOU NEED BINARY FORM OF SYMBOL AND 
DATA STATES INCLUDING SYMBOL OUTPUTS 
ENTER 2 IF YOU NEED DECIMAL FORM OF SYMBOL AND 
DATA STATES INCLUDING SYMBOL OUTPUTS 
ENTER 3 IF YOU WANT TO TURN BACK MAIN PROGRAM. 


1 
PLEASE ENTER THE SELECTION CORRECTLY 
ENTER 1 IF YOU NEED BINARY FORM OF SYMBOL AND 
DATA SIATES INCLUDING SYMBOL OUTPUTS 
ENTER 2 IF YOU NEED DECIMAL FORM OF SYMBOL AND 
DATA STATES INCLUDING SYMBOL OUTPUTS 
ENTER 3 IF YOU WANT TO TURN BACK MAIN PROGRAM. 


ANALYSES RESULTS WAS SENT TO THE FILE FTOQ8FOO1 
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~ 


SOURCE CODE 


THIS PROGRAM SIMULATES A M-DPSK COMMUNICATION SYSTEM 
PROGRAM RUNS TOTATLY INTERACTIVE THEN SENDS TO OUTPUTS FILES; 
Mber FrOGrOO] ,FILE RCDATA,FILE XMDATA,FILE R,FILE ERR. 

PRESSE PtLES ARE USED FOR INFORMATION TRANSFER BETWEEN 
SUBROUTINES BUT ARE NOT REQUIRED FOR PROGRAM EXECUTION. 

20U CAN SIMULATE &=-DPSK AND 16-DPSK BY USING THIS PROGRAM 
FOR 32-DPSK OR 64-DPSK SYSTEMS YOU HAVE TO CHANGE DIMENSIONS 
AND FORMAT STATEMENTS ACCORDINGLY,ALSO YOU HAVE TO EXPEND 
CHARACTER VARIABLE SET IN THE SUBROUTINE XMTR. 


OBERABBON SPECIFICATIONS....... ; 


THE SYMBOL STATES REQUESTEd DEPEND ON YOUR SYSTEM AND CAN 
Beall SEQUBNGE. MAIN THING IS; YOU HAVE TO SELECT YOUR REFERENCE 
Seweleneteem BY SONE STEP INCREMENTS, ALSO SAME SPECIFICATION 
eto ls POR YOUR ANGELES. PROGRAM ONLY EACEPIS THE INTEGER EQUIVALENTS 
OF YOUR STATE OR ANGLE CODES. 
eee eeu” CCC CCCECCECCCCCCCCCCCCECCCCECCECCECCCCCCEECCECECCCCCCCECCCCCCCC 


VARIABLE DECLERATION 


AAAAANAARANANAANANAANANANNMYA 


INTEGER N,L,ORDER,TM(64,64),PS(64) , TRDATA(256) ,M,KK,K, THETA(64) 
INTEGER MAT(64) , KL, RMAT(64,64) ,TR(64) ,TH(64,64) ,RCVMAT (64,64) ,NN 
INTEGER 90U OU, S U,LD,IS,NT(1,8),KM(1,8) ,KN(1,83 

INTEGER RCV(64,64) ,T(256) ,TD(256) ,NEXTST (256) , PRESTA(256) , FL 
INTEGER A,B,Y,KD,KY(1,8) ,KH(1,8) ,KJ(1,8) ,KB(1,8) 


CHARACTER*1 DATA(64,64) 
CHARACTER*16 SET 


é RK CALLING SUBROUITINES ARK 
CALL INFO (ORDER,PS,THETA,NN,N,FL) 
IF (FL.EQ.1)THEN 
CALL XMTR (ORDER, FL) 
CALL XMSTA (THETA,TM,PS,ORDER) 
CALL XMDAT(TM, ORDER, THETA ,NEXTST,PS,T,TD) 
CALL CONVER (N,NEXTST,NN,ORDER, FL) 
WRITE(6,1) 
1 FORMAT(T9,'’ ANALYSES RESULTS WAS SENT TO THE FILE FTO8FOO1') 


Meoe 
CALL RCVR(ORDER) 
CALL RCVSTA (PS,ORDER,THETA,RCVMAT) 
CALL RCVDAT(TR,PS,ORDER,NEXTST,RCVMAT ) 
CALL CONVER (N,NEXTST,NN,ORDER,FL) 
Se eer ore EO PORDER, TE TM tNeeret 1D) 
R ' 
2 FORMAT(T9,' ANALYSES RESULTS WAS SENT TO THE FILE FTO8FOO1',/, 
*T7,'ERROR ANALYSES RESULTS WAS SENT TO THE FILE ERR') 
ENDIF 
SLOP 
END 


SUBROUTINE INFO(ORDER,PS,THETA,NN,N,FL) 


THIS SUBROUTINE PROVIDES INFORMATIONS ABOUT OUR SYSTEM 


INTEGER N,L,ORDER,TM(64,64) ,PS(64) , THETA(64) , FL 
CHARACTER*1 DATA(64, 64) 


AQAA ee) 


WRITE (6,*) ' PLEASE SELECT ONE OF THEM! 
WRITE (6,*) ! ” IF YOU HAVE TRANSMITTER'| 
WRITE (6.*) '(2) IF YOU HAVE RECEIVER! 
READ ( A) FL 

WRITE (8. ' PLEASE SELECT ONE OF THEM! 
WRITE (6,*) '(1) IF YOUR SYSTEM 8-DPSK' 


89 


ANAAAA A 


AAAAM 


20 
uh 


10 
33 


81 


88 


WRITE 6 zx) '(2) IF YOUR SYSTEM 16-DPSK' 
READ (6,*) L 

IF(L.EO.1) THEN 

ORDER=8 

NN=6 3 

N=6 


ENDIF 

TE (lL EOe2Z oe ine 
ORDER=16 

N=8 
NN=255 
ENDIF 
WRITE 


‘ENTER THE INTEGER EQUIVALENTS OF SYMBOL STATES' 
aan THE ORDER OF COUNTER CLOCK WISE' 
R 


'ENTER THE INTEGER EQUIVALENTS OF THE ' 
‘PHASE ANGLES OF YOUR CODING SYSTEM ' 
ER 


iL 
,13,'*ALPHA= ?') 

READ(6,*) THETA(L) 

CONTINUE 

RETURN 


END 
SUBROUTINE XMTR(ORDER, FL) 


THIS SUBROUTINE PROVIDES GENERAL TRANSMITER MATRIX FOR 
M-DPSK SYSTEM. 
ITS OUTPUT CAN BE SEEN ON EITHER FILE FTO8F001 OR FILE 
XMDATA. 
INTEGER N,L,ORDER,FL 
CHARACTER*16 SET 
SET(1:16)= 'ABCDEFGHIJKLMNOP | 
OPEN(UNIT=11,FILE='XMDATA! , STATUS='OLD' ) 
IF (FL.EQ.1) THEN 
WRITE (8,81 ORDER 
FORMAT(///,T15, 'TRANSMITTER'S ENCODER MATRIX ANALYSIS 
* FOR' 13, '-DPSK' 
*,/,60('-'),//,75, ' GENERAL TRANSMITTER MATRIX! ,/,T2,31('-')) 


WRITE (8,88) 

pORIRE()// "TAs, TRANSUETTER INFORMATION FOR THE GIVEN RECEIVER! 
zx /,72('='),//,T75,'GENERAL TRANSMITTER MATRIX',/,T2,31('-')) 
ENDIF 


DO 3 N=1,ORDER 


IF (N.EQ. 1 THEN 

WRITE ,*) SET(N:ORDER) 

WRITE (11,*) SET(N:ORDER) 

ELSE 

WRITE ) SET(N:ORDER) ,SET(1:N-1) 
WRITE (11,*) SET(N:ORDER) ,SET(1:N-1) 
ENDIF 

CONTINUE 

CLOSE (11) 

RETURN 

END 


SUBROUTINE XMSTA(THETA,TM,PS,ORDER) 
THIS SUBROUTINE CONVERTS GENERAL XMTR MATRIX IN TO OUR SPECIAL 


AMTR SYMBOL MATRIA. 


GENERAL XMTR MATRIX IS READ FROM FILE XMDATA. 
INTEGER TM(64,64) ,PS(64) , THETA(64) ,ORDER 


90 


ANANA OY 


OCIA ©) 


og 


10 


43 


20 


83 
a 


77 


93 
ee 


~] OO 


CHARACTER*1 DATA(64,64) 
OPEN(UNIT=11,FILE='XMDATA! ,STATUS ='OLD') 
IF (ORDER EQ.8)THEN 
wees, 25) (THETA(T) ,I=1, ORDER) 
'TRANSMITTER MATRIX ',/,31('- a 
ota rs / DATA STATE (ROWS) 
pel, 19)" ( DArACE Ny N=158))1=1,8) 
FORMAT (1X, 8A1) 
ELSE 
WRITE (8 43) (THETA(T) , I=1 , ORDER) 


FORMAT (///,T 'TRANSMITTER MATRIX | Heo '), ar 
*'SYMBOL STATE ea. / DATA STATE (ROWS ) 
2 neh rt” im \ 
Bp (11,20) ((DATA(L, N),N=1,16),L=1,16) 
FORMAT (1X,16A1) 
ENDIF 


DO 6 M= 0,ORDER 

TE(T CHAR (DATA(K ND). .EQ. ICHAR(DATA(1, ORDER-M) ) ) THEN 
TH(K,N)=PS (ORDER-M 

ENDIF 

CONTINUE 

CONTINUE 

IF (ORDER.EQ.8)THEN 

WRITE (TM(K,N),N=1, ORDER) 

PS(K) ( TM(K,N) ,N=1, ORDER) 





ORDER) 
) ,N=1, ORDER) 


CONTINUE 
RETURN 
END 


SUBROUTINE XMDAT(TM,ORDER,THETA,NEXTST,PS,T,TD) 


THIS SUBROUTINE FINDS THE NEXT SYMBOL STATE/S BY DETERMINING 
foe UA cRIX CONSTRUCTED BY SUB. AMSTA 


INTEGER TM(64,64), ce 4), REE IW RET Oee NN 
INTEGER ORDER,T(256),TD(256), NEXTST (256), BPRESlA 2563 
DO 5 I=1, ORDER 

T(PS(I 1 

TD(THETA(I) )=I 

CONTINUE 

KK=0 

DO 7 K=0,ORDER-1 

DO 8 7 O (CRBER-1 

KK=KK+ 

NEXTST(RR)=TH(T(K), TD(M) ) 

CONTINUE 

CONTINUE 

RETURN 

END 


SUBROUTINE RCVR (ORDER) 


THIS SUBROUTINE PREPARES GENERAL RECEIVER MATRIX AND 
SEBS OUTPUTS TO THE "FILE RCDATA". 


INTEGER ORDER,MAT(64),KL,RMAT(64,64) 


OPEN (UNIT=2, FILE='RCDATA' , STATUS='OLD' ) 
WRITE (8,81)ORDER 


a 


81 


t4 


66 


77 
>> 


10 
i 


ue 
16 


88 


ANAAAAN A 


og 


£99 


vl 


eZ 


ig 
83 


x FOR',13,'-DPSK',/, 

SoOt =) / eS, ' GENERAL RECEIVER MATREX',/,T2,31('=-')) 
DO 3 I=0,0 DER-1 

MAT (I+1)=I 

CONTINUE 


DO 44 K=1,ORDER 

RMAT (1 ,K)=MAT(K) 
CONTIN 

DO 5 ee 2, ORDER 

DO 66 M=1,L- 

RMAT(L,M) “MAT (ORDER- L+M+1) 
CONTINUE 

DO 77 N=L,ORDER 

RMAT (L,N)=MAT (N-L+1) 
CONTINUE 
CONTINUE 

DO 88 NK=1,ORDER 

IF (ORDER.EQ.8)THEN 

WRITE (2. 10 (RHAT (NK, KL) Le) Oanee 
WRITE (8,15) (RMAT(NK,KL) ,KL=1,ORDER 
OE eetTs 

FORMAT (8X, 8(12)) 

ELSE 


WRITE er ET SU 
WRITE (8,16) (RMAT(NK,KL) ,KL=1,ORDER 
oa ey 

FORMAT (8X, 16(12)) 

ENDIF 

CONTINUE 

CLOSE(2) 

RETURN 

END 


SUBROUTINE RCVSTA (PS,ORDER,THETA,RCVMAT) 


THIS SUBROUTINE PREPARES SPECIFIC RECEIVER MATRIX UNDER 
DEFINED STATES.IT READS GENERAL RECEIVER MATRIX FORM FILE 
RCDATA WHICH HAD BEEN CREATED BY THE SUBROUTINE RCVDAT. 

OUTPUT ARE SENT TO THE FILE R. 


INTEGER TH(64,64) ,ORDER, PS(64) , THETA(64) ,RCVMAT(64,64) 
OPEN (ONTT= =2,FILE='RCDATA' , STATUS='OLD' ) 

OPEN(UNIT=20, FILE='R' -STATUS='! OLD! ) 

IF (ORDER.EQ.8)T 
WRITE (8. Sa) testeyat L=1,ORDER), (THETA (7 ) ,I=1,ORDER) 
FORMAT(///,1T5,'RECEIVER MATRIX 31('- af 

*| SYMBOL STATE (COLUMNS) / DATA eae (Rows)! , //,4X,8(I3)) 


ELS 
WRITE(S, 199 tre thy ,L=1,ORDER) 

FORMAT (///,T5,'RECEIVER MATRIX IS',/,31('- 3 ee 

*x'SYMBOL STATE (COLUMNS) / DATA STATE (ROWS)! 

* Ld AK, 16(I3) y 


DO. i T=1 

IF (ORDER. 20.8) 

READ (2, i Ht ON), N=1,ORDER) 
FORMAT (812 

ELSE 

READ (2, Be) (TH(I,N),N=1,ORDER) 
FORMAT (1612) 

ENDIF 

DO 2 K=1,ORDER 

Ra /K)=THETA(TH(I,K)+1) 


CONTIN 

IF (CRDERNE m: .8) THEN 

jaa 20, . (RCVMAT(I,N) ,N=1,ORDER) 
FORMAT 

WRITE a 83), i (RCVMAT(I,N) ,N=1,ORDER) 
FORMAT (I3 8(135) 


22 


ELSE 


Wh Te Sey PS(I),(RCVMAT(I,N) ,N=1,ORDER) 
93 Pomme (13,'|',16(I3)) 
WRITE(20,14) (RCVMAT(I,N) ,N=1,ORDER) 


14 FORMAT(1612) 


ENDIF 
1 CONTINUE 
CLOSE(20) 
RETURN 
2 END 
: SUBROUTINE RCVDAT(TR,PS,ORDER ,NEXTST,RCVMAT ) 
C THIS SUBROUTINE READS SPECIFIC RECEIVER MATRIX WHICH HAD 
€ BEEN CREATED BY SUBROUTINE RCVSTA,THEN SEARCHS FOR THE RCVR 
E OUTPUT. 
e 
é 
INTEGER RCVMAT(64,64) ,PS(64) ,KK,ORDER 
INTEGER TR(64) ,NEXTST(256) 
OPEN (UNIT=20, FILE='R! , STATUS='OLD' ) 
DO 23 LL=1,ORDER 
IF (ORDER.EO.8)THEN 
READ eg LL RM) /ie2 ORDER) 
24 FORMAT Baa 
& 
READ eee (LL, Ki) ,RM=1, ORDER) 
94  FORMAT(16(12)) 


ENDIF 
23 CONTINUE 

BO 2 I=1,ORDER 
M=PS (I) 
TR(M)=I 

Z CONTINUE 
KK=0 
BO 3 H=0,0RDER-1 
DO 4 K=0,ORDER-1 
KK=KK+1 
NEXTST(KK)=RCVMAT(TR(M) , TR(K) ) 

~ CONTINUE 

5 CONTINUE 


RETURN 
END 
E 
E 
; SUBROUTINE CONVER(N,NEXTST,NN,ORDER, FL) 
S 
INTEGER N,S,QU,SU,LLL,LLK(50) ,MMM,NN,FL,MM,M,K(9,1024) 
: INTEGER MW,NT(256,4) ,KZ,ORDER,NEXTST(256) ,ST,QQU,KD 


Sue Weiewe (6,*) ‘PLEASE ENTER THE SELECTION CORRECTLY ' 


WRITE (6,*%) ' ENTER 1 IF YOU NEED BINARY FORM OF SYMBOL AND | 
WRITE (6,*) 'DATA STATES INCLUDING SYMBOL OUTPUTS! 

WRITE (6,*) ' ENTER 2 IF YOU NEED DECIMAL FORM OF SYMBOL AND ' 
WRITE (6,%*) 'DATA STATES INCLUDING SYMBOL OUTPUTS' 

WRITE (6.*) ' ENTER 3 IF YOU WANT TO TURN BACK MAIN PROGRAM. | 
READ (6,*) S 

IF (S.EQ.1) GOTO 100 

Ip (SenO.2) GOTO 200 

IF (S.EO.3) GOTO 400 


IF (S.NE.1.0R.S.NE.2.0R.S.NE.3) GOTO 300 
200 IF (FL.EQ.1) THEN 
WRITE(8,6 
6 FORMAT (//,8X, ‘DATA STATE',5X,'SYMBOL STATE! ,5X, 'SYMBOL' ) 
ELSE 
WRITE(8,21) 
21 FORMAT (//,8X,'DATA STATE! ,5X,'SYMBOL STATE',5X,'SYMBOL' ) 


oS 


ANAAN A 


~J 00 


100 
47 


1a 


33 
45 


22 
45 
55 
46 


cS 
400 


DO 7 KZ=0,ORDER-1 

DO 8 M=0,ORDER-1 

KK=KK+1 

WRITE(8,9) KZ,M, pr 
FORMAT (10X,13,15X,13,15X, 13 
CONTINUE 

CONTINUE 

GOTO 300 

IF (FL.EQ.1) THEN 


FORMAT (//,7 a STATE! ,10X,'SYMBOL STATE! ,10X, 
*! SYMBOL, fe acre )) 


WRITE(8, oe 
FORMAT (//  SINBOL STATE! ,10X,'DATA STATE! ,10X, 
*'! SYMBOL | Le a3 


ENDIF 
ENTE a eo aytaew 
KD=3 


ELSE 
KD=4 

ENDIF 

DO 44 L=1,NN+1 

QOQU= SNEXTST(L) 

DO 33 SR=KD, 

NT(L, sk) =H05 (00, 2) 
QQU=09U/2 
CONTINUE 
CONTINUE 

DO 23 MW=0,NN 
QU=MW 

DO 22 LLL=N,1,-1 
K(MW, LLL)=HOD (QU, 2) 
QU=9U/2 
CONTINUE 

IF (ORDER.EQ. <8) THEN 


WRITE(8,45) MW, (K(MW, a ,LK=1,N), (NT(MW+1,ST),ST=1 
FORMAT (13, 2%,3(2%,12,2K), 3K, 32k 12 e) oe ee 
WRITE(8,55) MW, (K(MW, be ,LK=1,N), (NT (MW+1,ST),ST=1 a 
FORMAT (13,2X,4(2X,12,1X) ,3X,4(2X,12,1X) ,3X,4(2X,12,1X)) 

ENDIF 

WRITE (8. 46) 
FORMAT (72('=-')) 
CONTINUE 

GOTO 300 
RETURN 

END 


SUBROUTINE ERROR(FL,THETA,RCVMAT,PS,ORDER,TR,TM,T,NEXTST,TD) 
THIS SUBROUTINE LOCATES UNDETECTED ERROR COMBINATIONS. 


INTEGER QOQU, sO sou, ,PS(64),ORDER, RO, \oaees) Me a ote TM(64,64) 
INTEGER TD(256) ,T(256),A,B,Y,KD,KN(1,8),Kd(1,8),KB(1,8) ,KY(1,8) 
INTEGER FL,IS ee RACE, '8),NEXTST(256) , THETA(64) , KH(1,8) 


OPEN(UNIT=15,FILE='ERDATA' ,STATUS='OLD' ) 
WRITE (6,%) | IF YOU WANT TO HAVE UNDETECTABLE ERROR' 
WRITE x ‘ANALYSES ENTER "1" ,OTHERWISE ENTER "2"! 


LD 
aes (ORDER. EQ.8)THEN 


ELSE 
KD=4 
ENDIF 

IF( ED. ae x 
WRITE fee ; 
WRITE (6,* 


OPEN (UNT T=! - , FILE= 'ERR ,STATUS=' OLD! ) 


THEN 
WHAT KIND OF ANALYSES OUTPUT DO YOU WANT' 
(1) IN BINARY FORM! 
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iS 


eZ 


30 


Zee 


444 


So 


45 


85 


6 


35 


10 


6 A oe (2) IN DECIMAL FORM! 
READ AD a 





= 

1,72('#!),f/ IK, ‘CORRECT RECEPTION COMBINATIONS IN PRESE 
ROR',//, mot 41) 

T 


5, INPUT! ,T17,'PREVIOUS' ,T30,'PRESENT!' ,T43,'PREVIOUS', 
ENT! ,T68, 'OUTPUT',/,T15, ‘TRANSMISSION! ,T28, ‘TRANSMISSION’ 
* 242, ‘RECEPTION! ,T54, ‘RECEPTION! ./,79('-')) 


E 

Gomer l1 

ENDIF 

FL=2 

CALL XMTR (ORDER,FL) 

CALL XMSTA (THETA,TM,PS,ORDER) 
CALL XMDAT(TM,ORDER, THETA,NEXTST,PS,T,TD) 
DO 1 I=0,ORDER-1 

DO 2 K=0 ORDER~1 

A=TM(T(I) ,TD(K)) 

B= acumaT (ER (2). ee) 

eS ae 


DO 333 JJ=KD,1,-1 
NT(1, 3a) =10 (6u, 2) 
QQU=00U/2 

CONTINUE 

QU=K 

DO 222 LM=KD,1,-1 
ele /LM)=MOD(QU, 2) 


CONT NUE 

SQU=I 

DO 444 LN=KD,1,-1 
KN(1, ae =HOD(SQU, 2) 


oagussou/ E 
SQU=B 
DO 555 LN=KD, 
eh Le LN)= MoD (SOU, 2) 
SOU=SOU/2 
CONTINUE 
IF (ORDER.E 


.8) THEN 
ae, 45) KAI 
3) 


0) BU e SRE 
K= 


K Nin), 21-1, KD) 
De eee 
) 


* (1 IJ=1,KD 
FORMAT iToges (6s) Tac. BV es cicra 


RR Le CEC ie 


) (NT (1, 1T) IT=1,KD) 
15), TJ j= 
T 412), T65 ,4(I2)) 


Pisces loeloo,.loj.7),13) 


, ORD 
Ts (1 EQ. _f.AND.J.EQ.A)THEN 
Goro 1 

ZLSE 

= eee TAAL) TR(J)) 

Ti 

a BP Od) ty THEN 


SQU= 2 
DO 666 LN=KD,1,-1 
95 


666 


Tel, 


888 


48 


hele 


49 


949 


20 


100 


46 


Pre 


KH(1, By se a 2) 
SQU=SQU/ 
CONTINUE 
SQU=J 
DO 777 LN=KD,1,-1 
Koi. bw) =H0D (SOU, 2) 
SQU= =SOU We 
CONTINUE 
SQU=Y 
DO 888 LN=KD 
Kv, bw) =HoD (S6U, Zs 
SQU=SOU/2 
CONTINUE 
IF(ORDER.EO. 
WRITE(1,48 
Aer IT 
FORMAT (T4 


) 
H 


x OO 


RR) REEL, /KD), (KJ(1,KK) , KK=1, KD) 
T54,3(1I3),T65 3(73)) 
I ) 


iO~ = 
WO Ww “-_ 
| 


KK), KK=1, KD KJ(1,KK) ,KK=1,KD) 


> 
a 
ti 
J 
1 
~~ 
t+ 
i 
Hm" 1.0 to~ 
| 
ae | 


N~ 
x 


) ,KK=1,KD), (KJ(1,KK) , KK=1,KD) 


4,4(I2),T65 wea) 
KK) ,RK=1, KD), (KJ(1,KK) , KK=1, KD) 


LO Ar iad 
a. 


oN" OND” WO 
ri~ sR * WY 
hOM 
SA] 


Wand P20) eben y 
FORMAT (T48,13, £53. Eo, f7i ts) 
ENDIF 

ENDIF 
ENDIF 

CONTINUE 
CONTINUE 

WRITE (1,46) 
FORMAT (79('-')) 
CONTINUE 
CONTINUE 
CLOSE(1) 

Glesn (15) 
RETURN 

END 
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APPENDIX B 
COMPUTER GENERATED ENCODER AND DECODER OUTPUTS 


1 8-DPSK TRANSMITTER 


TRANSMITTER'S ENCODER MATRIX TABLE FOR 8-DPSK 


ABCDEFGH 
BCDEFGHA 
CDEFGHAB 
DEFGHABC 
EFGHABCD 
FGHABCDE 
GHABCDEF 
HABCDEFG 


TRANSMITTER MATRIX 


/ tap Ss) 2 toe 7 5 4 
i fee oles 2.0 64°66 
5 Hele soe 90) 4 6 7 
1 moe (2a 64/6 7 5 
3 Sec 084 6 7 5 1 
2 Paeomecwo F 5 1 3 
0 teeoeswoe / 5 61h U3 
fee olhUmflCUCUCOlUlUCUC dL ll CO” 
awe 7 5 132 0 4 
SYMBOL STATE DATA STATE SYMBOL 
0 0 0 0 0 0 0 0 0 ©) 
i 0 0 0 0 0 i I, 0 0 
Z 0 0 0 0 ft 0 i 1 1 
3 0 0 0 0 i 1 1 1 ©) 
4 0 0 0 1 0 0) 0 1 0 
2 0 0 0 t 0 d 0 t ih 
6 0 0 0 1 1 0 1 0 iE 
7 0 0 0 1 IL 1 0 0 1 
8 0 0 f 0 0 0 0 0 t 
9 ) 0 1 0 0 i 0 1 iE 
EG 0 0 I 0 i 0 ©) 0 0 
iL 0 0 1 0 i i 0 il 0 
EZ 0 0 1 1 0 0 1 0 i 


13 


0 G 


14 | 


0 


uf 


0 


1 


0 


Ea 


0 


18 


0 


iL 


20 


0 i 


0 


0 iL 


1 


0 1 0 


21 


iL 1 


i 


i 0 


1, 


0 i 0 


22 


0 1 


1 


1 1 


1 


1 0 


0 


23 


24 


i 0 


0 


0 1 


0 


1 1 


0 


25 


0 0 


di: 


iE 0 


0 


0 1 i 


26 


0 0 


0 


Za 


28 


1 


1 


Zo 


0 


0 


30 


i ib 


1 


1 Z 


1 


1 1 


0 


on 


0 


32 


1 0 


iL 


0 1 


0 


0 0 


if 


oS 


0 i. 


i 


1 0 


0 


q 0 0 


34 


BS 


0 0 


0 


0 0 


Ih 


1 0 0 


36 


i 0 


0 


0 L 


1 


0 0 


il 


a 


0 1 


0 


1 0 


1 


Ih 0 0 


38 


1 Z 


0 


1 1 


1 


V 0 0 


oo 


0 i 


i: 


0 0 


0 


0 1 


1 


40 


0 iE 


0 


0 1 


0 


0 i 


iL 


41 


0 


42 


i 


43 


1 


44 


45 


46 


47 


98 


48 


49 


50 


at 


Bic 


= 


ei 


56 


of 


58 


i 


50 


60 


61 


99 


8-DPSK RECEIVER 


RECEIVER'S DECODER MATRIX TABEE SHOR Bicep eon 


GENERAL RECEIVER MATRIX 


mOMNSTMONTAO 
OWNYFONRTAOE” 
MIHONAOMW 
TONTAOMr-WWM 
MNTAOM WM tH 
NADM™WOM HOH 
AODMUOMTON 
Orm~wuondtdrannre 


RECEIVER MATRIX IS 


SYMBOL STATE (COLUMNS) / SYMBOL STATE (ROWS) 


oS 20 4°36 


OFTM™ONMAO 


TNMmONMOTAOH 


MMmUONMNAOTM 


mUONMOAOTMHE 


ONMNAGDTMOS™ wo 


NO AOFTMO™UON 


MOnAOTFM™ONCH 


AOFM ONMr4 


OMmMAMNODSWO 


DATA STATE 


SYMBOL STATE 


SYMBOL STATE 


10 


ted 


wz 


3 


14 


ic 


100 


7 
i 


18 


0 


0 


19 


Ze 


25 


26 


re 


Z9 


0 


30 


gi 


a2 


38 


34 


36 


37 


38 


33 


40 


41 


0 


42 


43 


2¢ 


45 


46 


48 


4&9 


50 
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0 


0 


1 


TRANSMITTER INFORMATION FOR THE GIVEN RECEIVER 
5) 


3 26 7 


1 


GENERAL TRANSMITTER MATRIX 
0 


ABCDEFGH 
BCDEFGHA 
CDEFGHAB 
DEFGHABC 


EFGHABCD 
TRANSMITTER MATRIX 


GHABCDEF 
HABCDEFG 


Sye 
FGHABCDE 


54 


58 


> 


Sy 


62 


63 


60 


61 


SYMBOL STATE (COLUMNS) / DATA STATE (ROWS) 


OMmMNAMNOS 
POMMAMNO 
ODOM MTAMN 
NOTOrMAM 
MONOTORrMH 
raAMNO DHOr wm 
MAMNOAHAOr 


MmMTAMNODSWO 


MmUMNAMNODTWO 
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3s 16-DPSK TRANSMITTER 


TRANSMITTER'S ENCODER MATRIX TABLE FOR 8-DPSK 


ABCDEFGHIJKLMNOP 
BCDEFGHIJKLMNOPA 
CDEFGHIJKLMNOPAB 
DEFGHIJKLMNOPABC 
Er GHIJKLMNOPABCD 
r GHIJKLMNOPABCDE 
GHIJKLUNOPABCDEF 
HIJKLMNOPABCDEFG 
ITJKLMNOPABCDEFGH 
JKLUNOPABCDEFGHI 
KLMUNOPABCDEFGHIJ 
LMNOPABCDEFGHIJK 
MNOPABCDEFGHIJKL 
NOPABCDEFGHIJKLM 
OPASCDEFGHIJKLMN 
PABCDEFGHIJKLMNO 





TRANSMITTER MATRIX 


Weems © 7 > 4€ 1213 15 141011 9 8 


Plomts es tt 3 1 5 7 6 4 0 210 812 14 
mares 1 3 1 5S 7 6 4€ QO 210 8 12 14 15 
9 aoe oe >)? 6©66lUC CUCU 10 L2meol > 1s 
epee so lee > 6 f 6064 OO U622:10~=6812141513 9 
3 See > 4 One 100 dB Liz 14 15 13 9 11 
FE Meo ee ome 6 62 TO: 668 cdlzZ 614.1513 9 11 3 
2 ome fee eee 1006 68lLdi2 141513 911 3 1 
7 feo ewe 2 ro) 66612 614 15 13 8 1ll 3 1 5 
6 eee © 2 10 iets ls Soelt 3 Lt 5 7 
meeemrooeeetoe Ss 12 14 15 13 912 3 1 5 7 6 
0 Wee lOermert2 14915513 911 3 15 7 6 4 
Zz Pereeroeteetse 15 13 F911 3 1 5 7 6 4 OQ 
Metemecrizeta 15 13 911 3 1 5 7? 6 4 OO 2 
8 Smee tee torls 911 3 1 5 F? 6 4€ O 210 
oe tees Weis 911 3s 1 5 7? 6 4€ O 210 8 
Mave to is Fil 3 15 F 6 4 0 210 8 12 
a2meOL STATE DATA STATE SYMBOL 
0 0 0 0 0 0 0 0 0 0 0 0 0 
: 0 0 0 0 0 0 0 i, 0 0 i 0 
2 0 0 0 0 0 0 it 0 il 0 0 0 
3 0 0 0 0 0 0 1 1 1 0 1 0 
4 0 0 0 0 0 1 0 0 1 i 0 a 
5 0 0 0 0 0 I 0 x a a 1 1 
6 0 0 0 0 0 i 1 0 if 1 0 0 
iL 0 0 0 0 0 I I i Jl a 1 0 
8 0 0 0 0 2 0 0 0 0 I 0 0 


ip ett, mu a 


0 


10 


gk 


TZ 


oh 


14 


16 


7 


ES 


0 


19 


20 


0 


ZT 


0 


0 


22 


24 


26 


ay 


0 


28 


0 


29 


30 


Sel 


32 


35 


0 


34 


35 


36 


37 


0 


39 


40 


41 


42 


43 


104 


<j} 
<j' 


45 


46 


48 


4&9 


opi 


a2 


= 


55 


a0 


a 


58 


0 


a 


60 


61 


0 


62 


0 


63 


0 


0 


64 


66 


67 


1 


68 


0 


1 


0 


69 


0 


0 


0 


0 


70 


0 


0 


ve. 


0 


0 


72 


0 


13 


0 


74 


78 


105 


ie 


80 


81 


82 


83 


0 0 


0 


84 


86 


on 


88 


89 


90 


1 


a1 


92 


32 


94 


a0 


0 


96 


37 


96 


SF 


0 


100 
101 


102 


103 


104 


105 


106 


107 


108 


0 


16s 


Tro 


EIU Ik 


a2 


1 


Tes 


106 


5 


0 


i 


Tye? 


118 


11 


120 


Te 1 


Wz 2 


lize 


124 


125 


1 


1Ze6 


0 


17 


126 


0 


129 
130 


1 


ise 


136 


134 


135 


1 


136 


nS 7 


138 


130 


140 


141 


142 


143 


1 


144 


1 


145 


146 


147 


0 


148 


107 


149 


150 


| 


Ts2 


155 


154 


1 


135 


iSe 


iy 


158 


139 


160 


161 


ez 


163 


164 


ie 


166 


167 


168 


Le? 


170 


Ge 


it 


igzZ 


I 


is 


174 


175 


176 


i 


Tg 


Te 


Oe 


180 


Tel 


182 


183 


108 


184 


1 


185 


186 


1387 


188 


189 


190 


191 


PZ 


i 


194 


P25 


i 


196 


7 


198 


O9 


i 


200 


1 


201 


i 


202 


a 


203 


iL 


204 


1 


205 


i 


206 


207 


208 


209 


1 


210 


a 


i 


Ze 


1 


Z13 


i: 


214 


1 


a5 


1 


eaage@e@g@e@aq@aaae@a SBaeee@ataa@ae Bae @ eae S&F Bw ae este See 8B Bega at at ata at eee a SF SB eB et FB wat Ss Bees w_—sw st as SS Ss a = 


216 


I 


Za 


1 


218 


0 


0 
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0 


1 


ZA 
220 


0 


222 


0 


0 


22k 


223 


0 


1 


224 
225 


226 


Z 


2el 


223 


220 


250 


251 


1 


2a2 


233 


0 


0 


0 


0 


i 


234 


1 


235 


236 


2eu 


236 


0 


0 


0 


0 


1 


0 


0 


240 


0 


ca0 


241 


242 


243 


iE 


244 
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246 


247 


248 
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250 
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1 
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0 
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0 


1 
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4, 16-DPSK RECEIVER 


RECEIVER'S DECODER MATRIX TABLE FOR 16-DPSK 


OF PO GIP Un 


LOTTI Zeta , 
3101112131415 
9101112131415 

I1TOLT Zi siatS 
9101112131415 
6, 910 Zl si a henOn 
8 9ICl Ti 2 isso 
So 9G ligisilatS 


OrNWPROO~] 


em ho Ww & O10s~)] OC 
NG »s&10.~) 0 
WAH 305~) 0 


RECEIVER MATRIX IS 


15.13 9-1) 3 Pees 7 ee OO oe ee 
eS 0 1 3 2 6°77 5 42 AS ole Cee 
Nes 8 0 Ll 3 2-6 7 -5 42 128 ee 
g 9 8 O 1] -.3 2 6 2735) 482s oe 
li} 11°°9 8 Q° 1 3 260 64 eZee ee 
3; 10 11 9 8 OQ 1 3° 2,06 7 35 4a eee 
1} 14 10 ll 9 8 O 1/3 2°96 Hoes 
9| 15 14 10 11 9 §& OL 93) 2 5G peo eee 
7} 13.15 14 10 21 69 «8 00 ae 
6; 12.13 15 14 10 11-9 ¢ 01 eo eee 
4 41213 15 14 10 11 39° 8 GC) te Zoe 
0 9 41213 15 14 10 11) (9 6) Clase 
Z 7 5 41213 15 14 10) 119) 6 One 
10 6 7 § 4 12 13 15 14 165 Se yO ee ee 
8 2 6 7? 5 412 13 15 110 
LZ 3 2 6 7 5 412 13 15314 Croce ee 
i4h 1 3 2 6 7 § 4 12°13 15514 2001) oe 
SYMBOL STATE SYMBOL STATE DATA STATE 
0 0) 0 0) 0 0 0 0 () 0 0 0 0 
1 0 0 0 0) 0 0 0 1 i 1 1 0 
Z Q 0 0 0) 0 0 z 0 0 0) 0 1 
S 0 0 0 0 0 0 1 1 i I, a FE 
4 0 0 0 ) 0 I. 0 0 1 0 0 0 
5 0 0 0 0 0 1 0 1 1 Q 1 0 
6 0 0 0 0) 0 i 1 0 1 0 0 1 
d 0 0 0 0 0 1 it il i 0 if 1 
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S3 
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ABCDEFGHIJKLMNOP 
BCDEFGHIJKLMNOPA 
CDEFGHI JKLMNOPAB 
DEFGHIJKLMNOPABC 
EFGHI JKLMNOPABCD 
FGHI JKLMNOPABCDE 
GHIJKLMNOPABCDEF 
HIJKLMNOPABCDEFG 
IJKLMNOPABCDEFGH 
JKLMNOPABCDEFGHI 
KLMNOPABCDEFGHIJ 
LMNOPABCDEFGHIJK 
MNOPABCDEFGHIJKL 
NOPABCDEFGHIJKLM 
OPABCDEFGHIJKLMN 
PABCDEFGHIJKLMNO 


TRANSMITTER MATRIA 


15; 15 13 Sell 3 1 5 Fero 4€ 0 210 78 124 
13; 13 911.3 1 S92 6 42 GS Zete sete aS 
9 911 3 1 5° 7 GmeteeOo 2 105s eee 
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APPENDIX C 
COMPUTER GENERATED OUTPUTS FOR DBDR’S ERROR PATTERN 
ANALYSES 
The program outputs given below are a partial sample of the main outputs. 
Since undetectable error patterns are similar, in the next sections, 1/8 period of the 
pattern have been included. Actually 526 combinations of error patterns are detected 
by the program for the 8-DPSK case. For the 16-DPSK case there are 4096 


undetectable error pattern combinations. 
1. 8-DPSK RECEIVER 
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